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Abstract 

In these lectures the present status of the so-called standard cosmological model, 
based on the hot Big Bang theory and the inflationary paradigm is reviewed. 
Special emphasis is given to the origin of the cosmological perturbations we 
see today under the form of the cosmic microwave background anisotropics 
and the large scale structure and to the dark matter and dark energy puzzles. 



1 Introduction 

The evolution of the universe is determined to a large extent by the same microphysics laws of physics 
that govern high-energy physics phenomena. Hence, any progress in particle physics has a large impact 
on the cosmological model(s) and, conversely, any new step taken towards the understanding of the 
past, present and future of our universe might provide a hint of high-energy physics beyond the one we 
currently know. This is the reason why these lectures are entitled Particle Cosmology. If the reader 
takes only one lesson home from them it is that particle physics and cosmology are nowadays intimately 
connected. 

There are fundamental questions we are on the edge of answering: what is the origin of our 
universe? Why is the universe so homogeneous and isotropic on large scales? What are the origins of 
dark matter and dark energy? What is the fate of our universe? While these lectures will certainly not 
be able to give definite answers to them, we shall try to provide the students with some tools they might 
find useful in order to solve these overwhelming mysteries themselves. 

These lectures will contain a short review of the standard Big Bang model; a rather long discus- 
sion of the inflation paradigm with particular emphasis on the possibility that the cosmological seeds 
originated from a period of primordial acceleration; the physics of the Cosmic Microwave Background 
(CMB) anisotropics, and a discussion of the dark matter and dark energy puzzles. 

Since these lectures were delivered at a school, we shall not provide an exhaustive list of references 
to original material, but refer to several basic cosmology books and reviews where students can find the 
references to the original material |l]-[8j. 



2 Basics of the Big Bang model 

We know two basic facts about our local universe (the universe we may observe). First, it is homogeneous 
and isotropic on sufficiently large cosmological scales [2]. Once this experimental evidence is accepted, 
one can promote it to a principle, dubbed "the cosmological principle". Secondly, it expands. The next 
question would then be: how can we describe such a universe? 

The standard cosmology is based upon the maximally spatially symmetric Friedmann-Robertson- 
Walker (FRW) line element 
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where a(t) is the cosmic-scale factor, i? curv = a{t)\k\~ 1 / 2 is the curvature radius, and k = —1,0, 1 is 
the curvature signature. All three models are without boundary: the positively curved model is finite and 
curves back on itself; the negatively curved and flat models are infinite in extent. The Robertson-Walker 
metric embodies the observed isotropy and homogeneity of the universe. It is interesting to note that 



this form of the line element was originally introduced for the sake of mathematical simplicity; we now 
know that it is well justified at early times or today on large scales 10 Mpc), at least within our visible 
patch. 

The coordinates, r, 8, and (f>, are referred to as co-moving coordinates: A particle at rest in these 
coordinates remains at rest, i.e., constant r, 9, and <f). A freely moving particle eventually comes to rest 
in these coordinates, as its momentum is redshifted by the expansion, p oc a . Motion with respect to 
the co-moving coordinates (or cosmic rest frame) is referred to as peculiar velocity; unless supported by 
the inhomogeneous distribution of matter, peculiar velocities decay away as a -1 . Thus the measurement 
of peculiar velocities, which is not easy as it requires independent measures of both the distance and 
velocity of an object, can be used to probe the distribution of mass in the universe. 

Physical separations between freely moving particles scale as a(t); or said another way the physi- 
cal separation between two points is simply a(t) times the coordinate separation. The momenta of freely 
propagating particles decrease, or redshift, as a(t) , and thus the wavelength of a photon stretches as 
a(t), which is the origin of the cosmological redshift. The redshift suffered by a photon emitted from a 
distant galaxy 1 + z = ao/a(t); that is, a galaxy whose light is redshifted by 1 + z, emitted that light 
when the universe was a factor of (1 + z)~ l smaller. When the light from the most distant quasar yet 
seen (z = 4.9) was emitted, the universe was a factor of almost six smaller; when CMB photons last 
scattered, the universe was about 1100 times smaller. 



2.1 Friedmann equations 

The evolution of the scale factor a(t) is governed by Einstein equations 

R/J.U — 2 R9(iv — G»v = 8ttG , T MJy (2) 

where (p, v = 0, ■ ■ ■ 3) is the Riemann tensor and R is the Ricci scalar constructed via the metric 
|TJ) [2], and T^ u is the energy-momentum tensor. G = rripf is the Newton constant. Under the hypothesis 
of homogeneity and isotropy, we can always write the energy-momentum tensor under the form T^ v = 
diag (p, P, P, P) where p is the energy density of the system and P its pressure. They are functions of 
time. 

The evolution of the cosmic-scale factor is governed by the Friedmann equation 
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where p is the total energy density of the universe, matter, radiation, vacuum energy, and so on. 

Differentiating wrt to time both members of Eq. Q and using the the mass conservation equation 

p + 3H{p + P) = Q, (4) 

we find the equation for the acceleration of the scale factor 

a 6 

Combining Eqs. Q and ([5]) we find 

H = -4itG {p + P). (6) 
The evolution of the energy density of the universe is governed by 

d(pa 3 ) = -Pd (a 3 ) ; (7) 

which is the first law of thermodynamics for a fluid in the expanding universe. (In the case that the stress 
energy of the universe is comprised of several, non-interacting components, this relation applies to each 
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separately; e.g., to the matter and radiation separately today.) For P = p/3, ultra-relativistic matter, 
p oc a and a ~ £2 ; for P = 0, very nonrelativistic matter, p oc a and a ~ i 3 ; and for P = — p, 
vacuum energy, p = const. If the rhs of the Friedmann equation is dominated by a fluid with equation of 
state P = wp, it follows that p oc or 3(1+u,) and a oc t 2/ ' i{ - l+w \ 

We can use the Friedmann equation to relate the curvature of the universe to the energy density 
and expansion rate: 

0-1 = -^; « = («) 

a -n Pcrit 

and the critical density today p crit = 3H 2 /8irG = 1.88/i 2 gcm" ~ 3 ~ 1.05 x 10 4 eV cm . There is 
a one-to-one correspondence between SI and the spatial curvature of the universe: positively curved, 
Slo > 1; negatively curved, Slo < 1; and flat (Oo = 1). Further, the fate of the universe is determined by 
the curvature: model universes with k < expand forever, while those with k > necessarily recollapse. 
The curvature radius of the universe is related to the Hubble radius and S7 by 

TT-l 

Rcuiv = |n_i|i/2 ■ (9) 

In physical terms, the curvature radius sets the scale for the size of spatial separations where the effects 
of curved space become pronounced. And in the case of the positively curved model it is just the radius 
of the 3-sphere. 

The energy content of the universe consists of matter and radiation (today, photons and neutrinos). 
Since the photon temperature is accurately known, To = 2.73 ± 0.01 K, the fraction of critical density 
contributed by radiation is also accurately known: fi/j/i 2 = 4.2 x 10~ 5 , where h = 0.72 ± 0.07 is 
the present Hubble rate in units of 100 km s _1 Mpc -1 [9]. The remaining content of the universe 
is another matter. Rapid progress has been made recently toward the measurement of cosmological 
parameters [ [TO] ]. Over the past years the basic features of our universe have been determined. The 
universe is spatially flat; accelerating; comprised of one third dark matter and two thirds a new form of 
dark energy. The measurements of the cosmic microwave background anisotropics at different angular 
scales performed by the WMAP Collaboration Q have recently significantly increased the case for 
accelerated expansion in the early universe (the inflationary paradigm) and at the current epoch (dark 
energy dominance), especially when combined with data on high-redshift supernovae (SN1) and large- 
scale structure (LSS) (TO). The CMB+LSS+SN1 data give (9j 

f2 = 1.00+o;o3 ) 

meaning that the present universe is spatially flat (or at least very close to being flat). Restricting to 
S7q = 1, the dark matter density is given by (9j 

O h 2 — r\ 1 1+0.0034 
"DM" — u - li -0.059 ' 

and a baryon density 

n B = 0.045 ± 0.0015, 

while the Big Bang nucleosynthesis estimate is ilsh 2 = 0.019 ± 0.002. Substantial dark (unclustered) 
energy is inferred: 



Ode « 0.72 ±0.015. 

What is most relevant for us is that this universe was apparently born from a burst of rapid expansion, 
inflation, during which quantum noise was stretched to astrophysical size seeding cosmic structure. This 
is exactly the phenomenon we want to address in part of these lectures. 
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2.2 The early, radiation-dominated universe 

In any case, at present, matter outweighs radiation by a wide margin. However, since the energy density 
in matter decreases as a~ 3 , and that in radiation as a~ 4 (the extra factor due to the redshifting of the 
energy of relativistic particles), at early times the universe was radiation dominated — indeed the calcu- 
lations of primordial nucleosynthesis provide excellent evidence for this. Denoting the epoch of matter 
and radiation equality by subscript 'EQ,' and using To = 2.73 K, it follows that 

a EQ = 4.18 x 1(T 5 (fWr 1 ; T EQ = 5.62(f} /j 2 ) eV ; (10) 

t EQ =4.17x W 10 (n h 2 )- 2 s. (11) 

At early times the expansion rate and age of the universe were determined by the temperature of the 
universe and the number of relativistic degrees of freedom: 

Prad = <?*(T) — ; H ~ 1.67 5 yV/mpi; (12) 

=^aoct 1/2 ; t ~ 2.42 x 10"V 1/2 (^/GeV)- 2 s; (13) 

where g* (T) counts the number of ultra-relativistic degrees of freedom (a* the sum of the internal degrees 
of freedom of particle species much less massive than the temperature) and mpi = G -1 / 2 = 1.22 x 
10 19 GeV is the Planck mass. For example, at the epoch of nucleosynthesis, g* = 10.75 assuming three, 
light (^C MeV) neutrino species; taking into account all the species in the Standard Model, g* = 106.75 
at temperatures much greater than 300 GeV. 

A quantity of importance related to g* is the entropy density in relativistic particles, 

T 45 y ' 

and the entropy per co-moving volume, 

S oc a 3 s oc 5*a 3 T 3 . 

By a wide margin most of the entropy in the universe exists in the radiation bath. The entropy density is 
proportional to the number density of relativistic particles. At present, the relativistic particle species are 
the photons and neutrinos, and the entropy density is a factor of 7.04 times the photon-number density: 

n 7 = 413 cm -3 and s = 2905 cm~ 3 . 

In thermal equilibrium — which provides a good description of most of the history of the universe — 
the entropy per co-moving volume S remains constant. This fact is very useful. First, it implies that the 
temperature and scale factor are related by 

roc <7*~ 1/3 a-\ (14) 

which for g* = const leads to the familiar T oc a -1 . 

Second, it provides a way of quantifying the net baryon number (or any other particle number) per 
co-moving volume: 

N B = R\ B = ^ ~ (4 - 7) x 10 -11 . (15) 

The baryon number of the universe tells us two things: (1) the entropy per particle in the universe is 
extremely high, about 10 10 or so compared to about 10~ 2 in the Sun and a few in the core of a newly 
formed neutron star. (2) The asymmetry between matter and antimatter is very small, about 10~ 10 , since 
at early times quarks and antiquarks were roughly as abundant as photons. One of the great successes of 
particle cosmology is baryogenesis, the idea that B, C, and CP violating interactions occurring out-of- 
equilibrium early on allow the universe to develop a net baryon number of this magnitude. Finally, the 
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constancy of the entropy per co-moving volume allows us to characterize the size of co-moving volume 
corresponding to our present Hubble volume in a very physical way: by the entropy it contains, 

Su = fH^s * 10 90 . (16) 

The standard cosmology is tested back to times as early as about 0.01 s; it is only natural to ask how 
far back one can sensibly extrapolate. Since the fundamental particles of Nature are point-like quarks and 
leptons whose interactions are perturbatively weak at energies much greater than 1 GeV, one can imagine 
extrapolating as far back as the epoch where general relativity becomes suspect, i.e., where quantum 
gravitational effects are likely to be important: the Planck epoch, t ~ 10 _43 s and T ~ 10 19 GeV. Of 
course, at present, our firm understanding of the elementary particles and their interactions only extends 
to energies of the order of 100 GeV, which corresponds to a time of the order of 10~ n s or so. We can 
be relatively certain that at a temperature of 100-200 MeV (t ~ 10~ 5 s) there was a transition (likely a 
second-order phase transition) from quark/gluon plasma to very hot hadronic matter, and that some kind 
of phase transition associated with the symmetry breakdown of the electroweak theory took place at a 
temperature of the order of 300 GeV (f ~ 10~ n s). 



2.3 The concept of particle horizon 

In spite of the fact that the universe was vanishingly small at early times, the rapid expansion precluded 
causal contact from being established throughout. Photons travel on null paths characterized by dr = 
dt/a(t); the physical distance that a photon could have travelled since the bang until time t, the distance 
to the particle horizon, is 



Ra(t) ~ a(t) I W 



H 

n 



(1-n) (1-n 



for a(t) oaf 1 , n < 1. 



Using the conformal time dr = dt/a, the particle horizon becomes 



R H (t) = a(r) / dr, 

J TQ 



(17) 



(18) 



where to indicates the conformal time corresponding to t = 0. Note, in the standard cosmology the 
distance to the horizon is finite, and up to numerical factors, equal to the age of the universe or the 
Hubble radius, H -1 . For this reason, we shall use horizon and Hubble radius interchangeably 1 . 

Note also that a physical length scale A is within the horizon if A < Rh ~ H -1 . Since we can 
identify the length scale A with its wavenumber k, A = 2ira/k, we shall have the following rule 



_k_ 

aH 

_k_ 

all 



< 1 SCALE A OUTSIDE THE HORIZON 



> l SCALE A WITHIN THE HORIZON 



1 As we shall see, in inflationary models the horizon and Hubble radius are not roughly equal as the horizon distance grows 
exponentially relative to the Hubble radius; in fact, at the end of inflation they differ by e N , where iV is the number of e-folds of 
inflation. However, we shall slip and use "horizon" and "Hubble radius" interchangeably, though we shall always mean Hubble 
radius. 
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3 The shortcomings of the standard Big Bang theory 

By now the shortcomings of standard cosmology are well appreciated: the horizon or large-scale smooth- 
ness problem; the small-scale inhomogeneity problem (origin of density perturbations); and the flatness 
or oldness problem, we shall briefly review only the horizon problem here here. 



3.1 The horizon problem 

According to standard cosmology, photons decoupled from the rest of the components (electrons and 
baryons) at a temperature of the order of 0.3 eV. This corresponds to the so-called surface of 'last scat- 
tering' at a redshift of about 1100 and an age of about 180000 (f) fr 2 ) _1/2 yr- From the epoch of last 
scattering onwards, photons free-stream and reach us basically untouched. Detecting primordial pho- 
tons is therefore equivalent to take a picture of the universe when the latter was about 300 000 years 
old. The spectrum of the cosmic background radiation (CBR) is consistent with that of a black body at 
temperature 2.73 K over more than three decades in wavelength. 

The most accurate measurement of the temperature and spectrum is that by the WMAP5 instru- 
ment on the COBE satellite which determined its temperature to be 2.726 ± 0.01 K (9J. The length 
corresponding to our present Hubble radius (which is approximately the radius of our observable uni- 
verse) at the time of last scatteringwas 

A/KtLs) = R H {h) (—) = R H (t ) ( 



On the other hand, during the matter-dominated period, the Hubble length decreased with a different law 

H 2 oc pm °c a~ 3 oc T 3 . 



At last-scattering 



Hll = R H (to) (^A ' <R H (to). 



The length corresponding to our present Hubble radius was much larger that the horizon at that time. 
This can be by shown comparing the volumes corresponding to these two scales 



hi! " v^* 



2 ss 10 6 . (19) 



There were ~ 10 6 casually disconnected regions within the volume that now corresponds to our horizon! 
It is difficult to come up with a process other than an early hot and dense phase in the history of the 
universe that would lead to a precise black body for a bath of photons which were causally disconnected 
the last time they interacted with the surrounding plasma. 

The horizon problem is well represented by Fig. [T] where the solid line indicates the horizon scale 
and the dashed line any generic physical length scale A. Suppose, indeed, that A indicates the distance 



between two photons we detect today. From Eq. (19) we discover that at the time of emission (last- 



scattering) the two photons could not talk to each other, the dashed line is above the solid line. There is 
another aspect of the horizon problem which is related to the problem of initial conditions for the cos- 
mological perturbations. We have every indication that the universe at early times, say t <C 300 000 yr, 
was very homogeneous; however, today inhomogeneity (or structure) is ubiquitous: stars (5p/ p ~ 10 30 ), 
galaxies (5p/p ~ 10 5 ), clusters of galaxies (5p/p ~ 10 — 10 3 ), superclusters, or "clusters of clusters" 
(5p/ p ~ 1), voids (5p/ p ~ —1), great walls, and so on. For some twenty-five years standard cosmology 
has provided a general framework for understanding this picture. Once the universe becomes matter 
dominated (around 1000 yr after the bang) primeval density inhomogeneities (Sp/p ~ 10 -5 ) are ampli- 
fied by gravity and grow into the structure we see today (2j. The existence of density inhomogeneities 
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log a 



Fig. 1: The horizon scale (solid line) and a physical scale A (dashed line) as function of the scale factor a 
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Fig. 2: The CMBR anisotropy as function of £ (from Ref. |9|) 



has another important consequence: fluctuations in the temperature of the CMB radiation of a similar 
amplitude. The temperature difference measured between two points separated by a large angle (> 1°) 
arises due to a very simple physical effect: the difference in the gravitational potential between the two 
points on the last scatteringsurface, which in turn is related to the density perturbation, determines the 
temperature anisotropy on the angular scale subtended by that length scale, 



ST\ 
~Tj, 



S£ 

P J A 



(20) 



where the scale A ~ lOO/i 1 Mpc(0/deg) subtends an angle 9 on the last-scattering surface. This is 
known as the Sachs-Wolfe effect |TT}[T2j. We shall come back to this piece of physics. The temperature 
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anisotropy is commonly expanded in spherical harmonics 



— (x ,t , n) = ^ a e,m(x )Y em (n), (21) 

em 

where xq and tq are our position and the preset time, respectively, n is the direction of observation, £'s 
are the different multipoles and 2 

( a em a *e' m >} = fie,e'$m,m'Ce, (22) 

where the deltas are due to the fact that the process that created the anisotropy is statistically isotropic. 
The Cg's are the so-called CMB power spectrum. For homogeneity and isotropy, the Cg's are neither a 
function of xq, nor of m. The two-point correlation function is related to the Cfs in the following way 

/5T(n) 5T(n') x ^ 



= E^E y ^( n ) y ^( n/ ) = ^E(^ +1 ) c '^(' u = n - n/ ) w 

em I 
where we have used the addition theorem for the spherical harmonics, and Pi is the Legendre polynom 



of order t. In expression ( 23 1 the expectation value is an ensemble average. It can be regarded as an 
average over the possible observer positions, but not in general as an average over the single sky we 
observe, because of the cosmic variance 3 . WMAP5 data are given in Fig. [2] 

Let us now consider the last scatteringsurface. In co-moving coordinates the latter is 'far' from us 
a distance equal to 

f to dt r° 

/ - = / dT = (T -T LS ). (24) 

Jt hS a Jt LS 

A given co-moving scale A is therefore projected on the last scatteringsurface sky on an angular scale 

A (25) 



(to - tls 



where we have neglected tiny curvature effects. Consider now that the scale A is of the order of the co- 
moving sound horizon at the time of last-scattering, A ~ c s tls> where c s ~ 1/ v3 is the sound velocity 
at which photons propagate in the plasma at the last-scattering. This corresponds to an angle 



6~c s - — — — -~c s — , (26) 

(TO-TLSj TO 

where the last passage has been performed knowing that to 3> tls- Since the universe is matter- 
dominated from the time of last scatteringonwards, the scale factor has the following behaviour: a ~ 

t 2/3 

~ t . The angle #hor subtended by the sound horizon on the last-scattering surface then 
becomes ^ 

#HOR - c s ~ 1°, (27) 

where we have used Tls — 0.3 eV and Tq ~ 10~ 13 GeV. This corresponds to a multipole ^hor 



2 An alternative definition is Ci = (|ei£ m | 2 ) = 2I+1 X] 



2t+l i—im.= -l 



3 The usual hypothesis is that we observe a typical realization of the ensemble. This means that we expect the difference 
between the observed values \ai m \ 2 and the ensemble averages Ct to be of the order of the mean-square deviation of |e^ m j 2 
from Ct. The latter is called cosmic variance and, because we are dealing with a Gaussian distribution, it is equal to 2Ce for 
each multipole I. For a single i, averaging over the (21 + 1) values of m reduces the cosmic variance by a factor {21 + 1), but 
it remains a serious limitation for low multipoles. 
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Fig. 3: The behaviour of a generic scale A and the horizon scale H 1 in the standard inflationary model 



^HOR = ~ Q 200 . (28) 

#HOR 

From these estimates we conclude that two photons which on the last scatteringsurface were separated 
by an angle larger than #hor> corresponding to multipoles smaller than ^hor ~ 200, were not in causal 
contact. On the other hand, from Fig. [2] it is clear that small anisotropics, of the same order of mag- 
nitude 5T/T ~ 10 -5 are present at £ <C 200. We conclude that one of the striking features of the 
CMB fluctuations is that they appear to be non-causal. Photons at the last scatteringsurface which were 
causally disconnected have the same small anisotropics! The existence of particle horizons in the stan- 
dard cosmology precludes explaining the smoothness as a result of microphysical events: the horizon at 
decoupling, the last time one could imagine temperature fluctuations being smoothed by particle inter- 
actions, corresponds to an angular scale on the sky of about 1°, which precludes temperature variations 
on larger scales from being erased. To account for the small-scale lumpiness of the universe today, den- 
sity perturbations with horizon-crossing amplitudes of 10~ 5 on scales of 1 Mpc to 10 4 Mpc or so are 
required. 

As can be seen in Fig. [TJ in the standard cosmology the physical size of a perturbation, which 
grows as the scale factor, begins larger than the horizon and, relatively late in the history of the uni- 
verse, crosses inside the horizon. This precludes a causal microphysical explanation for the origin of the 
required density perturbations. 

From the considerations made so far, it appears that solving the horizon problem of the standard 
Big Bang theory requires that the universe go through a primordial period during which the physical 
scales A evolve faster than the horizon scale H . 

If there is period during which physical length scales grow faster than H , length scales A which 
are within the horizon today, A < H^ 1 (such as the distance between two detected photons) and were 
outside the horizon for some period, A > H^ 1 (for instance at the time of last scattering when the two 
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photons were emitted), had a chance to be within the horizon at some primordial epoch, A < H again, 
see Fig. 3. If this happens, the homogeneity and the isotropy of the CMB can easily be explained: pho- 
tons that we receive today and were emitted from the last scattering surface from causally disconnected 
regions have the same temperature because they had a chance to 'talk' to each other at some primordial 
stage of the evolution of the universe. 

The second condition can easily be expressed as a condition on the scale factor a. Since a given 
scale A scales like A ~ a and H~ l = a /a, we need to impose that there is a period during which 

Gft) = i>0 ' 

We can therefore introduce the following rigorous definition: an inflationary stage is a period of the 
universe during which the latter accelerates 



INFLATION <=^> a > 0. 



Comment : Let us stress that during such an accelerating phase the universe expands adiabatically. 
This means that during inflation one can exploit the usual FRW equations ([3]) and Q. It must be clear 
therefore that the non-adiabaticity condition is satisfied not during inflation, but during the phase transi- 
tion between the end of inflation and the beginning of the radiation-dominated phase. At this transition 
phase a large entropy is generated under the form of relativistic degrees of freedom: the Big Bang has 
taken place. 

4 The standard inflationary universe 

From the previous section we have learned that an accelerating stage during the primordial phases of the 
evolution of the universe might be able to solve the horizon problem. From Eq. (|5]) we learn that 

a > + 3P) < . 

An accelerating period is obtainable only if the overall pressure p of the universe is negative: p < —p/3. 
Neither a radiation-dominated phase nor a matter-dominated phase (for which p = p/3 and p = 0, 
respectively) satisfy such a condition. Let us postpone for the time being the problem of finding a 
'candidate' able to provide the condition P < —p/3. For sure, inflation is a phase of the history of 
the universe occurring before the era of nucleosynthesis (t « 1 s, T m 1 MeV) during which the light 
elements abundances were formed. This is because nucleosynthesis is the earliest epoch from which we 
have experimental data and they are in agreement with the predictions of the standard Big Bang theory. 
However, the thermal history of the universe before the epoch of nucleosynthesis is unknown. 

In order to study the properties of the period of inflation, we assume the extreme condition p = —p 
which considerably simplifies the analysis. A period of the universe during which P = — p is called the 
de Sitter stage. By inspecting Eqs. ([3]) and Q, we learn that during the de Sitter phase 

p = constant , 
Hi = constant , 

where we have indicated by Hj the value of the Hubble rate during inflation. Correspondingly, solving 
Eq. Q gives 

a = a i e Hl{ - t - ti) , (29) 
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where ti denotes the time at which inflation starts. Let us now see how such a period of exponential 
expansion takes care of the shortcomings of the standard Big Bang Theory 4 . 



4.1 Inflation and the horizon problem 

During the inflationary (de Sitter) epoch the horizon scale H^ 1 is constant. If inflation lasts long enough, 
all the physical scales that have left the horizon during the radiation-dominated or matter-dominated 
phase can re-enter the horizon in the past: this is because such scales are exponentially reduced. As we 
have seen in the previous section, this explains both the problem of the homogeneity of CMB and the 
initial condition problem of small cosmological perturbations. Once the physical length is within the 
horizon, microphysics can act, the universe can be made approximately homogeneous and the primeval 
inhomogeneities can be created. 

Let us see how long inflation must be sustained in order to solve the horizon problem. Let U and 
tf be, respectively, the time of beginning and end of inflation. We can define the corresponding number 
of e-foldings N 

N = \n[H I (t e -t i )}. (30) 

A necessary condition to solve the horizon problem is that the largest scale we observe today, the present 
horizon Hq 1 , was reduced during inflation to a value \h {U) smaller than the value of horizon length 
HJ 1 during inflation. This gives 

w«-<£)(£)-V$).-»£*r'. 

where we have neglected for simplicity the short period of matter-domination and we have called Tf the 
temperature at the end of inflation (to be indentified with the reheating temperature Trh at the beginning 
of the radiation-dominated phase after inflation, see later). We get 

^KI)-KIH 7+in (l)- 

Apart from the logarithmic dependence, we obtain N > 70. 



4.2 A prediction of inflation 

Since during inflation the Hubble rate is constant 

o 1 k 1 

On the other hand it is easy to show that to reproduce a value of (Oo — 1) of order of unity today, the 
initial value of (fi — 1) at the beginning of the radiation-dominated phase must be \Sl — 1| ~ 10~ 60 . 
Since we identify the beginning of the radiation-dominated phase with the beginning of inflation, we 
require 

\n-i\ t=tf ~w- m . 

During inflation 

l^-i| t=ti W/ 

Taking |fi — l| t=t . of order unity, it is enough to require that N rj 70. However, IF the period of inflation 
lasts longer than 70 e-foldings the present-day value of S7o will be equal to unity with great precision. 
One can say that a generic prediction of inflation is that 



4 Despite the fact that the growth of the scale factor is exponential and the expansion is superluminal, this is not in contra- 
diction with what is dictated by relativity. Indeed, it is the spacetime itself which is progating so fast and not a light signal in 
it. 
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INFLATION 



fin 



1. 



This statement, however, must be taken cum grano salis and properly specified. Inflation does not change 
the global geometric properties of the space-time. If the universe is open or closed, it will always remain 
flat or closed, independently from inflation. What inflation does is to magnify the radius of curvature 
i? curv defined in Eq. Q so that locally the universe is flat with a great precision. As we shall see, the 
current data on the CMB anisotropics confirm this prediction. 



4.3 Inflation and the inflaton 



In the previous subsections we have described the various advantages of having a period of accelerating 
phase. The latter required P < — p/3. Now, we would like to show that this condition can be attained by 
means of a simple scalar field. We shall call this field the inflaton (p. 

The action of the inflaton field reads 



d 4 x v '—g £ 



d 4 x- 



-g 



dpWt + v{4>) 



(32) 



where yj— g = a 6 for the FRW metric (lib. From the Euler-Lagrange equations 



If 1 



8(^£) 5(y/=g£) 
8d^(f) 5<j> 



0, 



we obtain 



+ m<t) 



+ V'(<j>) = 0. 



(33) 



(34) 



where V'{4>) = (dV(<fi)/d<j)). Note, in particular, the appearance of the friction term 3H(j): a scalar field 
rolling down its potential suffers a friction due to the expansion of the universe. 

We can write the energy momentum tensor of the scalar field 

= d^cj)d u (j) - g^y £ . 

The corresponding energy density and pressure density P^ are 

Too = P4> = y + V ^> + ~2^2~' 



Pn 



6a 2 



(35) 



(36) 



Note that, if the gradient term were dominant, we would obtain P^ 
We can now split the inflaton field in 

4>{t) = 0o(t) + <fy(x,i), 



'- , not enough to drive inflation. 



where <po is the 'classical' (infinite wavelength) field, that is the expectation value of the inflaton field on 
the initial isotropic and homogeneous state, while 5(f)(x, t) represents the quantum fluctuations around 
4>o- In this section, we shall be concerned only with the evolution of the classical field 0o- The next 
section will be devoted to the crucial issue of the evolution of quantum perturbations during inflation. 
This separation is justified by the fact that quantum fluctuations are much smaller than the classical value 
and therefore negligible when looking at the classical evolution. Not to be overwhelmed by the notation, 
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we shall indicate the classical value of the inflaton field by <fi from now on. The energy momentum tensor 
becomes 



Too = P<s> = y + V{4>) (37) 
T ii = P <i> = ^-V{ct>). (38) 

If 

v(<f>) » 4> 2 

we obtain the following condition 

P<t> a -PJ, . 

From this simple calculation, we realize that a scalar field whose energy is dominant in the universe and 
whose potential energy dominates over the kinetic term gives inflation. Inflation is driven by the vacuum 
energy of the inflaton field. 



4.4 Slow-roll conditions 

Let us now quantify better under which circumstances a scalar field may give rise to a period of inflation. 
The equation of motion of the field is 



+ 3H<j> + V'( 



0. 



(39) 



If we require that 4> 2 <C V (</>)> tne scalar field is slowly rolling down its potential. This is the reason why 
such a period is called slow-roll. We may also expect that since the potential is flat, <p is negligible as 
well. We shall assume that this is true and we shall quantify this condition soon. The FRW equation Q 
becomes 



(40) 



where we have assumed that the inflaton field dominates the energy density of the universe. The new 
equation of motion becomes 

3H4> = -V'{<f>) (41) 



which gives as a function of V'(4>). Using Eq. (41 1 slow-roll conditions then require 



F « V{4>) 



V 



< H 



and 



< 3H<]) => V" < H 2 . 
It is now useful to define the slow-roll parameters e and 77 in the following way 



V 

s 



H 

1 fV" 

~v 



4vrG 



1 



8vrG 



H 2 16ttG 
_ 1 V" 
~ 3^2' 



/\ 2 



H(j) 
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It might be useful to have the same parameters expressed in terms of conformal time 



S = 7] — e = 1 — 



The parameter e quantifies how much the Hubble rate H changes with time during inflation. Notice that, 
since 



H + H 2 = (1 - e) H 2 , 



inflation can be attained only if e < 1 : 



INFLATION 



e < 1. 



As soon as this condition fails, inflation ends. In general, slow-roll inflation is attained if e <C 1 and 
\rj\ <C 1. During inflation the slow -roll parameters e and 77 can be considered to be approximately 
constant since the potential V((f>) is very flat. 

Comment : In the following, we shall work at first-order perturbation in the slow -roll parameters, 
that is we shall take only the first power of them. Since, using their definition, it is easy to see that 
e, fj = O (e 2 , T] 2 ) , this amounts to saying that we shall treat the slow -roll parameters as constant in time. 

Within these approximations, it is easy to compute the number of e-foldings between the beginning 
and the end of inflation. If we indicate by 4>i and 4> / the values of the inflaton field at the beginning and 
at the end of inflation, respectively, we find that the total number of e-foldings is 

/**/ 

N = Hdt 

hi 

* Hf*' d -t 

Hi 4> 

~ -8vrG / — #. (42) 
We may also compute the number of e-foldings AN which are left to go to the end of inflation 



AN ~ 8vrG 



r<t>i 

hi 



4>AN y 



7#, 



(43) 



where 4>an is the value of the inflaton field when there are AN e-foldings to the end of inflation. 

1. Comment : According to the criterion given in Subsection 2.4, a given scale length A = a/k 
leaves the horizon when k = aH^ where is the the value of the Hubble rate at that time. One can 
easily compute the rate of change of H 2 as a function of k 



dlrvH 2 
din k 



dlnH 2 
dt 



dt 



din a j \ din k 



din a 



H 1 

2— x — x 1 
H H 



'IP 



-2e. 



(44) 
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2. Comment : Take a given physical scale A today which crossed the horizon scale during inflation. 
This happened when 



a {I) ^ - a (M) e_A " - 1171 



where AN\ indicates the number of e-foldings from the time the scale crossed the horizon during infla- 
tion and the end of inflation. This relation gives a way to determine the number of e-foldings to the end 
of inflation corresponding to a given scale 

/ A \ / I/ 1 / 4 \ / T 

A7V A -65 + ln — — - +21n — u — - - In ' f 



3000 Mpc J \ 10 14 GeV J V 1Ql ° GeV 

Scales relevant for the CMB anisotropics correspond to AiV ~60. 

Inflation ended when the potential energy associated with the inflaton field became smaller than 
the kinetic energy of the field. By that time, any pre-inflation entropy in the universe had been inflated 
away, and the energy of the universe was entirely in the form of coherent oscillations of the inflaton 
condensate around the minimum of its potential. The universe may be said to be frozen after the end of 
inflation. We know that somehow the low-entropy cold universe dominated by the energy of coherent 
motion of the <p field must be transformed into a high-entropy hot universe dominated by radiation. The 
process by which the energy of the inflaton field is transferred from the inflaton field to radiation has 
been dubbed reheating. In the theory of reheating, the simplest way to envisage this process is if the co- 
moving energy density in the zero mode of the inflaton decays into normal particles, which then scatter 
and thermalize to form a thermal background. It is usually assumed that the decay width of this process 
is the same as the decay width of a free inflaton field. 

Of particular interest is a quantity usually known as the reheat temperature, denoted as Tr# 5 . The 
reheat temperature is calculated by assuming an instantaneous conversion of the energy density in the 
inflaton field into radiation when the decay width of the inflaton energy, T^, is equal to H, the expansion 
rate of the universe. 

The reheat temperature is calculated quite easily. After inflation the inflaton field executes co- 
herent oscillations about the minimum of the potential. Averaged over several oscillations, the coherent 
oscillation energy density redshifts as matter: p^ oc a~ 3 , where a is the Robertson-Walker scale factor. 
If we denote as pj and ai the total inflaton energy density and the scale factor at the initiation of coherent 
oscillations, then the Hubble expansion rate as a function of a is 

H\a) = 8 ^J ) L- 2 ( a L) 3 . (45) 
3 mpi z V a J 

Equating H{a) and leads to an expression for aj/a. Now if we assume that all available coherent 
energy density is instantaneously converted into radiation at this value of aj/a, we can find the re- 
heat temperature by setting the coherent energy density, p^ = pj(ai/a) 3 , equal to the radiation energy 
density, pr = (tt 2 /30)g*T^ H , where g* is the effective number of relativistic degrees of freedom at 
temperature Trh- The result is 

/ 90 X 1 / 4 /200\ 1//4 
TrH= W^) V^=0.2[-) yr>^. (46) 

5 Inflation and the cosmological perturbations 

As we have seen in the previous section, the early universe was made very nearly uniform by a primordial 
inflationary stage. However, the important caveat in that statement is the word 'nearly'. Our current un- 
derstanding of the origin of structure in the universe is that it originated from small 'seed' perturbations, 



5 So far, we have indicated it by Tj. 
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which over time grew to become all of the structure we observe. Once the universe becomes matter dom- 
inated (around 1000 yrs after the bang) primeval density inhomogeneities (Sp/p ~ 10 -5 ) are amplified 
by gravity and grow into the structure we see today [4]. The fact that a fluid of self- gravitating particles is 
unstable to the growth of small inhomogeneities was first pointed out by Jeans and is known as the Jeans 
instability. Furthermore, the existence of these inhomogeneities was confirmed by the COBE discov- 
ery of CMB anisotropics; the temperature anisotropics detected almost certainly owe their existence to 
primeval density inhomogeneities, since, as we have seen, causality precludes microphysical processes 
from producing anisotropics on angular scales larger than about 1°, the angular size of the horizon at 
last-scattering. 

The growth of small matter inhomogeneities of wavelength smaller than the Hubble scale (A < 
H~ l ) is governed by a Newtonian equation: 

k 2 

<5 k + 2H5± + v 2 -J <5 k = 47rG/>M<5k, (47) 

where v 2 = dP/dpM is the square of the speed of sound and we have expanded the perturbation to the 
matter density in plane waves 

Wx,t) - 1 fd 3 k5 k (t)e- lk *. (48) 



Pm (2vr) 3 

Competition between the pressure term and the gravity term on the rhs of Eq. ( |47| ) determines whether or 
not pressure can counteract gravity: perturbations with wavenumber larger than the Jeans wavenumber, 
k 2 = ATiGa 2 p m /v 2 , are. Jeans stable and just oscillate; perturbations with smaller wavenumber are Jeans 
unstable and can grow. 

Let us discuss solutions to this equation under different circumstances. First, consider the Jeans 
problem, evolution of perturbations in a static fluid, i.e., H = 0. In this case Jeans unstable perturbations 
grow exponentially, 5^ oc exp(i / r) where r = 1/ y/AGitpM- Next, consider the growth of Jeans unstable 
perturbations in a matter-dominated universe, i.e., H 2 = 87rGp M /3and a oc i 2 / 3 . Because the expansion 
tends to pull particles away from one another, the growth is only power law, 5^ oc i 2 / 3 ; i.e., at the same 
rate as the scale factor. Finally, consider a radiation-dominated universe. In this case, the expansion is 
so rapid that matter perturbations grow very slowly, as In a in a radiation-dominated epoch. Therefore, 
perturbations may grow only in a matter-dominated period. Once a perturbation reaches an overdensity 
of order unity or larger it separates from the expansion, i.e., it becomes its own self-gravitating system 
and ceases to expand any further. In the process of virial relaxation, its size decreases by a factor of 
two — density increases by a factor of 8; thereafter, its density contrast grows as a 3 since the average 
matter density is decreasing as a -3 , though smaller scales could become Jeans unstable and collapse 
further to form smaller objects of higher density. 

In order for structure formation to occur via gravitational instability, there must have been small 
pre-existing fluctuations on physical length scales when they crossed the Hubble radius in the radiation- 
dominated and matter-dominated eras. In the standard Big Bang model these small perturbations have to 
be put in by hand, because it is impossible to produce fluctuations on any length scale while it is larger 
than the horizon. Since the goal of cosmology is to understand the universe on the basis of physical laws, 
this appeal to initial conditions is unsatisfactory. The challenge is therefore to give an explanation to the 
small seed perturbations which allow the gravitational growth of the matter perturbations. 

Our best guess for the origin of these perturbations is quantum fluctuations during an inflation- 
ary era in the early universe. Although originally introduced as a possible solution to the cosmological 
conundrums such as the horizon, flatness and entropy problems, by far the most useful property of infla- 
tion is that it generates spectra of both density perturbations and gravitational waves. These perturbations 
extend from extremely short scales to scales considerably in excess of the size of the observable universe. 

During inflation the scale factor grows quasi-exponentially, while the Hubble radius remains al- 
most constant. Consequently the wavelength of a quantum fluctuation — either in the scalar field whose 
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potential energy drives inflation or in the graviton field — soon exceeds the Hubble radius. The amplitude 
of the fluctuation therefore becomes 'frozen in'. This is quantum mechanics in action at macroscopic 
scales. 

According to quantum field theory, empty space is not entirely empty. It is filled with quantum 
fluctuations of all types of physical fields. The fluctuations can be regarded as waves of physical fields 
with all possible wavelenghts, moving in all possible directions. If the values of these fields, averaged 
over some macroscopicaily large time, vanish then the space filled with these fields seems to us empty 
and can be called the vacuum. 

In the exponentially expanding universe the vacuum structure is much more complicated. The 
wavelenghts of all vacuum fluctuations of the inflaton field (f> grow exponentially in the expanding uni- 
verse. When the wavelength of any particular fluctuation becomes greater than H -1 , this fluctuation 
stops propagating, and its amplitude freezes at some non-zero value 5<p because of the large friction term 
3H(p the equation of motion of the field (f>. The amplitude of this fluctuation then remains almost un- 
changed for a very long time, whereas its wavelength grows exponentially. Therefore, the appearance of 
such frozen fluctuation is equivalent to the appearance of a classical field 5<f) that does not vanish after 
having averaged over some macroscopic interval of time. Because the vacuum contains fluctuations of all 
possible wavelengths, inflation leads to the creation of more and more new perturbations of the classical 
field with wavelength larger than the horizon scale. 

Once inflation has ended, however, the Hubble radius increases faster than the scale factor, so the 
fluctuations eventually re-enter the Hubble radius during the radiation- or matter-dominated eras. The 
fluctuations that exit around 60 e-foldings or so before reheating re-enter with physical wavelengths in the 
range accessible to cosmological observations. These spectra provide a distinctive signature of inflation. 
They can be measured in a variety of different ways including the analysis of microwave background 
anisotropics. 

Quantum fluctuations of the inflaton field are generated during inflation. Since gravity talks to any 
component of the universe, small fluctuations of the inflaton field are intimately related to fluctuations 
of the space-time metric, giving rise to perturbations of the curvature 1Z (which will be defined in the 
following; the reader may loosely think of it as a gravitational potential). The wavelengths A of these 
perturbations grow exponentially and leave the horizon soon when A > Rh- On superhorizon scales, 
curvature fluctuations are frozen in and may be considered as classical. Finally, when the wavelength 
of these fluctuations re-enters the horizon, at some radiation- or matter-dominated epoch, the curvature 
(gravitational potential) perturbations of the space-time give rise to matter (and temperature) perturba- 
tions 5p via the Poisson equation. These fluctuations will then start growing, giving rise to the structures 
we observe today. 

In summary, these are the key ingredients for understanding the observed structures in the universe 
within the inflationary scenario: 

- Quantum fluctuations of the inflaton field are excited during inflation and stretched to cosmological 
scales. At the same time, being the inflaton fluctuations connected to the metric perturbations 
through Einstein's equations, ripples on the metric are also excited and stretched to cosmological 
scales. 

- Gravity acts a messenger since it communicates the small seed perturbations to photons and 
baryons once a given wavelength becomes smaller than the horizon scale after inflation. 

Let us now see how quantum fluctuations are generated during inflation, we shall proceed by 
steps. First, we shall consider the simplest problem of studying the quantum fluctuations of a generic 
scalar field during inflation: we shall learn how perturbations evolve as a function of time and compute 
their spectrum. Then — since a satisfactory description of the generation of quantum fluctuations has to 
take both the inflaton and the metric perturbations into account — we shall study the system composed 
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by quantum fluctuations of the inflaton field and quantum fluctuations of the metric. 



6 Quantum fluctuations of a generic massless scalar field during inflation 

Let us first see how the fluctuations of a generic scalar field x, which is not the inflaton field, behave 
during inflation. To warm up we first consider a de Sitter epoch during which the Hubble rate is constant. 

6.1 Quantum fluctuations of a generic massless scalar field during a de Sitter stage 

We assume this field to be massless. The massive case will be analysed in the next subsection. 
Expanding the scalar field x in Fourier modes 

we can write the equation for the fluctuations as 

k 2 

5xk + 3H5x k + ^5 Xk = 0. (49) 

or 



Let us study the qualitative behaviour of the solution to Eq. ( 49 1 



For wavelengths within the horizon, A <C H , the corresponding wave-number satisfies the 



relation k^> aH.ln this regime, we can neglect the friction term 3H <5Ak an d Eq. ( 49 1 reduces to 



k 2 

<%k + -2 <5*k = 0, (50) 

which is basically the equation of motion of an harmonic oscillator. Of course, the frequency 
term k 2 /a 2 depends upon time because the scale factor a grows exponentially. On the qualitative 
level, however, one expects that when the wavelength of the fluctuation is within the horizon, the 
fluctuation oscillates. 

For wavelengths above the horizon, A S> H , the corresponding wave-number satisfies the rela- 



tion k -C aH and the term k 2 /a 2 can be safely neglected. Equation (49 1 reduces to 

5xk + 3H6 Xk = 0, (51) 
which tells us that on superhorizon scales #Xk remains constant. 

We have therefore the following picture: take a given fluctuation whose initial wavelength A ~ 
a I k is within the horizon. The fluctuations oscillate till the wavelength becomes of the order of the 
horizon scale. When the wavelength crosses the horizon, the fluctuation ceases to oscillate and gets 
frozen in. 

Let us now study the evolution of the fluctuation in a more quantitative way. To do so, we perform 
the following redefinition 

, <fok 

o*k = 

a 

and we work in conformal time dr = dt/ a. For the time being, we solve the problem for a pure de Sitter 
expansion and we take the scale factor exponentially growing as a ~ e Ht ; the corresponding conformal 
factor reads (after choosing properly the integration constants) 

a(r) = --i- (r<0). 
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In the following we shall also solve the problem in the case of quasi de Sitter expansion. The beginning 
of inflation coincides with some initial time Tj <C 0. We find that Eq. d49b becomes 



5a'k + ( k 2 — — ) (Vr k 



0. 



(52) 



We obtain an equation which is very 'close' to the equation for a Klein-Gordon scalar field in flat space- 
time, the only difference being a negative time-dependent mass term —a" /a = — 2/r 2 . Equation (52 1 
can be obtained from an action of the type 



6Su 



dr 



K 2 



k 



a 
a 



5a 2 



(53) 



which is the canonical action for a simple harmonic oscillator with canonical commutation relations 

5 at 5ai - 5a^5at! = -i. 



Let us study the behaviour of this equation on subhorizon and superhorizon scales. Since 

_k_ 
~a~H 



-k t . 



on subhorizon scales k 2 ^> a" /a Equation (52) reduces to 

5a'k + k 2 <5cJk = . 

whose solution is a plane wave 



-ikr 



&7 k 



2 k 



(k » aH) . 



(54) 



We find again that fluctuations with wavelength within the horizon oscillate exactly like in flat space- 
time. This does not come as a surprise. In the ultraviolet regime, that is for wavelengths much smaller 
than the horizon scale, one expects that approximating the space-time as flat is a good approximation. 

On superhorizon scales, k 2 <C a" /a Equation (52i reduces to 



Say 



0, 



which is satisfied by 



Say, = B(k) a (k < aH) 



(55) 



where B{k) is a constant of integration. Roughly matching the (absolute values of the) solutions (54) 

de 

1 



and (55) at k = aH (—kr = 1), we can determine the (absolute value of the) constant B{k) 

H 



\B{k)\a 



\B(k)\ 



1 



'2k ' " " aV2k V2k 3 
Going back to the original variable 5xk, we obtain that the quantum fluctuation of the \ fi e ld on super- 
horizon scales is constant and approximately equal to 



I^Xkl - 



H 



(ON SUPERHORIZON SCALES) 



In fact we can do much better, since Eq. ( 52 ) has an exact solution: 

„— ikr 



5cr k 



1 + 



kr 



(56) 



This solution reproduces all that we have found by qualitative arguments in the two extreme regimes 
k <C aH and k ^> aH. We have performed the matching procedure to show that the latter can be very 
useful to determine the behaviour of the solution on superhorizon scales when the exact solution is not 
known. 
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6.2 The power spectrum 

Let us define now the power spectrum, a useful quantity to characterize the properties of the perturba- 
tions. For a generic quantity g(x, t), which can expanded in Fourier space as 

/r7 3 k 

the power spectrum can be defined as 

(Ol^kjO) = 5^ ( kl - k 2 ) ^P g (k), (57) 
where 1 ) is the vacuum quantum state of the system. This definition leads to the usual relation 

/dk 
-j-V g (k). (58) 

6.3 Quantum fluctuations of a generic scalar field in a quasi de Sitter stage 

So far, we have computed the time evolution and the spectrum of the quantum fluctuations of a generic 
scalar field x supposing that the scale factor evolves like in a pure de Sitter expansion, a(r) = —1/(Ht). 
However, during inflation the Hubble rate is not exactly constant, but changes with time as H = — e H 2 
(quasi de Sitter expansion). In this subsection, we shall solve for the perturbations in a quasi de Sitter 
expansion. Using the definition of the conformal time, one can show that the scale factor for small values 
of e becomes 

The fluctuation mass-squared mass term is 



M 2 (r) =m 2 x a 2 



a 
a 



where 

a" 



a 2 (^ + H ^j=a 2 (H + 2H 2 ) 
a 2 (2-e)H 2 - (2 " e) 



r 2 (l-ef 

~ ^(2 + 3e). (59) 

Armed with these results, we may compute the variance of the perturbations of the generic \ field 

d 3 k , . 2 



/d 6 k 
(2^)3 l<5xk| 



/ 



dk k 3 . ,2 
~k 2^ l<5xk| 

dk 

T V Sx (k), (60) 



which defines the power spectrum of the fluctuations of the scalar field \ 



V8 x {k) = ^ \5 X *\ 2 ■ (61) 
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Since we have seen that fluctuations are (nearly) frozen in on superhorizon scales, a way of characterizing 
the perturbations is to compute the spectrum on scales larger than the horizon. For a massive scalar field, 
we obtain 



V Sx (k) 



H Y ( k 



2vr / \aH 



where, taking m x /H 2 = 3r] x and expanding for small values of e and 77, 



(62) 



(63) 



We may also define the spectral index n$ x of the fluctuations as 



n 5x 



din Vxrh 

1 = - 7 = 3 - 2u x = 2n x - 2e. 
din k x ix 



The power spectrum of fluctuations of the scalar field x is therefore nearly flat, that is is nearly indepen- 
dent of the wavelength A = ir/k: the amplitude of the fluctuation on superhorizon scales does almost not 
depend upon the time at which the fluctuation crosses the horizon and becomes frozen in. The small tilt 
of the power spectrum arises from the fact that the scalar field x is massive and because during inflation 
the Hubble rate is not exactly constant, but nearly constant, where 'nearly' is quantified by the slow -roll 
parameters e. Adopting the traditional terminology, we may say that the spectrum of perturbations is 
blue if ns x > 1 (more power in the ultraviolet) and red if n$ x < 1 (more power in the infrared). The 
power spectrum of the perturbations of a generic scalar field x generated during a period of slow-roll 
inflation may be either blue or red. This depends upon the relative magnitude between 7] x and e. 

Comment : We might have computed the spectral index of the spectrum Vs x (k) by first solving the 
equation for the perturbations of the field x in a di Sitter stage, with H = constant and therefore e = 0, 
and then taking into account the time evolution of the Hubble rate introducing the subscript in whose 



time variation is determined by Eq. (44i. Correspondingly, is the value of the Hubble rate when a 
given wavelength ~ k^ 1 crosses the horizon (from that point on the fluctuation remains frozen in). The 
power spectrum in such an approach would read 



VsJk) 



H k \' 2 ( k ^ 3_2 "\ 



2vr J \aH 



(64) 



with 3 — 2v x ~ rj x . Using Eq. (44i, one finds 



n Sx 



dlnVss din Hi 
1 = 7 = —— * + 3 - 2v x = 2t] x - 2e 
dink M - u 



din k 



which reproduces our previous findings. 

Comment : Since on superhorizon scales 



<$Xk 



H 

fZk? 



k 



H 



1 + (Vx ~ e ) ln 



k 
a~H 



we discover that 



I^Xkl ^ \H(r) x - e) 5xk\ < \H 5xk\ , 
that is, on superhorizon scales the time variation of the perturbations can be safely neglected. 



(65) 



21 



7 Quantum fluctuations during inflation 

As we have mentioned in the previous section, the linear theory of the cosmological perturbations repre- 
sents a cornerstone of modern cosmology and is used to describe the formation and evolution of structures 
in the universe as well as the anisotropics of the CMB. The seeds for these inhomogeneities were gener- 
ated during inflation and stretched over astronomical scales because of the rapid superluminal expansion 
of the universe during the (quasi) de Sitter epoch. 

In the previous section we have already seen that pertubations of a generic scalar field \ WQ 
generated during a (quasi) de Sitter expansion. The inflaton field is a scalar field and, as such, we 
conclude that inflaton fluctuations will be generated as well. However, the inflaton is special from the 
point of view of perturbations. The reason is very simple. By assumption, the inflaton field dominates the 
energy density of the universe during inflation. Any perturbation in the inflaton field means a perturbation 
of the stress energy momentum tensor 



ST, 



A perturbation in the stress energy momentum tensor implies, through Einstein's equations of motion, a 
perturbation of the metric 



ST, 



fJ,V 



5R 



2<5 {g^vR) 



S-kGST, 



fJ,V 



Sg, 



in/- 



On the other hand, a pertubation of the metric induces a back-reaction on the evolution of the inflaton 
perturbation through the perturbed Klein-Gordon equation of the inflaton field 



flU 



5 £W + 



dV 







This logic chain makes us conclude that the perturbations of the inflaton field and of the metric are tightly 
coupled to each other and have to be studied together 



Sg^v ■ 



As we shall see shortly, this relation is stronger than one might think because of the issue of gauge 
invariance. 

Before launching ourselves into the problem of finding the evolution of the quantum perturbations 
of the inflaton field when they are coupled to gravity, let us give a heuristic explanation of why we expect 
that during inflation such fluctuations are indeed present. 

If we take Eq. ( [34] > and split the inflaton field as its classical value 4>q plus the quantum flucutation 
Scj), </>(x, t) = (f>o(t) + <5</>(x, t), the quantum perturbation 5(j) satisfies the equation of motion 



+ 3H, 



+ V 



0. 



(66) 



Differentiating Eq. ( 39 1 wrt time and taking H constant (de Sitter expansion) we find 



+ 3#^o + V" ■ 



0. 



(67) 



Let us consider for simplicity the limit k 2 /a 2 <C 1 and let us disregard the gradient term. Under this 
condition we see that 4>q and 5(f) solve the same equation. The solutions have therefore to be related to 
each other by a constant of proportionality which depends upon time 
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5(j) = —0o St(x). 



(68) 



This tells us that <p(x, t) will have the form 

(j>(x,t) = 0o (x,i- St(x)) . 

This equation indicates that the inflaton field does not acquire the same value at a given time t in 
all the space. On the contrary, when the inflaton field is rolling down its potential, it acquires different 
values from one spatial point x to the next. The inflaton field is not homogeneous and fluctuations are 
present. These fluctuations, in turn, will induce fluctuations in the metric. 

7.1 The metric fluctuations 

The mathematical tool to describe the linear evolution of the cosmological perturbations is obtained by 
perturbing at the first order the FRW metric gffl , see Eq. jlj 

9nu = 9$(t) + p^(x,t); gp,<g$. (69) 

The metric perturbations can be decomposed according to their spin with respect to a local rotation of 
the spatial coordinates on hypersurfaces of constant time. This leads to 

— scalar perturbations 

— vector perturbations 

— tensor perturbations 

Tensor perturbations or gravitational waves have spin 2 and are the true degrees of freedom of 
the gravitational fields in the sense that they can exist even in the vacuum. Vector perturbations are 
spin 1 modes arising from rotational velocity fields and are also called vorticity modes. Finally, scalar 
perturbations have spin 0. 

Let us do a simple exercise to count how many scalar degrees of freedom are present. Take a space- 
time of dimensions D = n + 1, of which n coordinates are spatial coordinates. The symmetric metric 
tensor g^ v has |(n + 2)(n + 1) degrees of freedom. We can perform (n + 1) coordinate transformations 
in order to eliminate (n+ 1) degrees of freedom, this leaves us with ^n(n+l) degrees of freedom. These 
\n(n+\) degrees of freedom contain scalar, vector and tensor modes. According to Helmholtz's theorem 
we can always decompose a vector m (i = 1, • • • , n) as Ui = div + V{, where v is a scalar (usually 
called potential flow) which is curl-free, vum = 0, and V{ is a real vector (usually called vorticity) 
which is divergence-free, V ■ v = 0. This means that the real vector (vorticity) modes are (n — 1). 
Furthermore, a generic traceless tensor Uij can always be decomposed as Ily = Hfj + 11^ + Hfj, where 

Hg = (-^p 2 + n, nYj = (-i/2k) (^n,- + %ILJ (^II; = 0) and ^Tlg = 0. This means that 

the true symmetric, traceless and transverse tensor degreees of freedom are i(n — 2)(n + 1). 

The number of scalar degrees of freedom is therefore 

^n(n + 1) - (n - 1) - ~(ra - 2)(n + 1) = 2, 

while the degrees of freedom of true vector modes are (n — 1) and the number of degrees of freedom of 
true tensor modes (gravitational waves) is |(n — 2)(n + 1). In four dimensions n = 3, meaning that one 
expects 2 scalar degrees of freedom, 2 vector degrees of freedom and 2 tensor degrees of freedom. As we 
shall see, to the 2 scalar degrees of freedom from the metric, one has to add another one, the inflaton field 
perturbation 5cf). However, since Einstein's equations will tell us that the two scalar degrees of freedom 
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from the metric are equal during inflation, we expect a total number of scalar degrees of freedom equal 
to 2. 

At the linear order, the scalar, vector, and tensor perturbations evolve independently (they decou- 
ple) and it is therefore possible to analyse them separately. Vector perturbations are not excited during 
inflation because there are no rotational velocity fields during the inflationary stage, we shall analyse the 
generation of tensor modes (gravitational waves) in the following. For the time being we want to focus 
on the scalar degrees of freedom of the metric. 

Considering only the scalar degrees of freedom of the perturbed metric, the most generic perturbed 
metric reads 

2 ( -1 - 2$ d t B \ 

9ilv ~ a V (1 - 2Y0<% + DijE J ' (/U) 

while the line-element can be written as 

ds 2 = a 2 {{-I -2<5>)cLt 2 + 2d i Bdrdx i + ((1 -2^)% + D^E) cteW). (71) 
Here = (didj - g^V 2 ). 



7.2 The issue of gauge invariance 

When studying the cosmological density perturbations, what we are interested in is following the evolu- 
tion of a space-time which is neither homogeneous nor isotropic. This is done by following the evolution 
of the differences between the actual space-time and a well understood reference space-time. So we shall 
consider small perturbations away from the homogeneous, isotropic space-time. 

The reference system in our case is the spatially flat Friedmann-Robertson-Waiker (FRW) space- 
time, with line element ds 2 = a 2 (r) {dr 2 — 5ijdx l dx^. Now, the key issue is that general relativity 
is a gauge theory where the gauge transformations are the generic coordinate transformations from one 
local reference frame to another. 

When we compute the perturbation of a given quantity, this is defined to be the difference between 
the value that this quantity assumes on the real physical space-time and the value it assumes on the 
unperturbed background. Nonetheless, to perform a comparison between these two values, it is necessary 
to compute them at the same space-time point. Since the two values live on two different geometries, 
it is necessary to specify a map which allows one to link univocally the same point on the two different 
space-times. This correspondence is called a gauge choice and changing the map means performing a 
gauge transformation. 

Fixing a gauge in general relativity implies choosing a coordinate system. A choice of coordinates 
defines a threading of space-time into lines (corresponding to fixed spatial coordinates x) and a slicing 
into hypersurfaces (corresponding to fixed time r). A choice of coordinates is called a gauge and there 
is no unique preferred gauge 



GAUGE CHOICE SLICING AND THREADING 



From a more formal point of view, operating an infinitesimal gauge transformation on the coordinates 

xP = x v + 5x^ (72) 
implies on a generic quantity Q a transformation on its perturbation 

SQ = 5Q + £ Sx Q (73) 
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where Qq is the value assumed by the quantity Q on the background and £$ x is the Lie-derivative of Q 
along the vector Sx^. 

Decomposing in the usual manner the vector 5x^ 

5x° = £°o M ) ; 

6x { = &P(af) + v^x^) ; 8^ = 0, (74) 

we can easily deduce the transformation law of a scalar quantity / (like the inflaton scalar field <fi and 
energy density p). Instead of applying the formal definition ( f73j ), we find the transformation law in 
an alternative (and more pedagogical) way. We first write Sf(x) = f(x) — fo(x), where fo(x) is the 
background value. Under a gauge transformation we have 6f(x^) = f(x^) — fo(x^). Since / is a scalar 
we can write f(xP) = f{x^) (the value of the scalar function in a given physical point is the same in 
all the coordinate system). On the other side, on the unperturbed background hypersurface /q = /q. We 
have therefore 



5f(x") 



f&) - /o(x?) 



/(^)-«5x'^(x)-/o(^) 



(75) 



from which we finally deduce, being /o = fo(x ), 



sf = sf-n 



For the spin-zero perturbations of the metric, we can proceed analogously. We use the following trick. 
Upon a coordinate transformation x M — > x^ = x^ + 5x^, the line element is left invariant, ds 2 = ds 2 . 



2 <x°) 



1 + $) (dx°) = a 2 (x°) (1 + $) (dx ) 2 . Since a 2 (x°) 



dx u + dx\ we obtain 1 + 2<I> = 1 + 2$ + 2 + 2^ u/ . We 



dx 



This implies, for instance, that a 

a 2 (x°) + 2a a' ^° and dx° = (l + £ 0/ 
now may introduce in detail some gauge-invariant quantities which play a major role in the computation 
of the density perturbations. In the following we shall be interested only in the coordinate transformations 
on constant time hypersurfaces and therefore gauge invariance will be equivalent to independence of the 
slicing. 



7.3 The co-moving curvature perturbation 

The intrinsic spatial curvature on hypersurfaces on constant conformal time r and for a flat universe is 
given by 

a z 

The quantity ip is usually referred to as the curvature perturbation. We have seen, however, that the 
curvature potential ip is not gauge invariant, but is defined only on a given slicing. Under a transformation 
on constant time hypersurfaces t — > t + 5t (change of the slicing) 
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tp ^ ip+ HSt. 

We now consider the co-moving slicing which is defined to be the slicing orthogonal to the worldlines of 
co-moving observers. The latter are are free-falling and the expansion defined by them is isotropic. In 
practice, what this means is that there is no flux of energy measured by these observers, that is Toi = 0. 
During inflation this means that these observers measure 5(j) com = since Tqi goes like <9j50(x, t)4>'(t). 

Since Sep — > 5(f> — 4>'5t for a transformation on constant time hypersurfaces, this means that 



5<f> — > 5<fi c 



b'5r = 



5t 



that is 5t = is the time-displacement needed to go from a generic slicing with generic 5(f) to the 
co-moving slicing where 54> CO m = 0. At the same time the curvature perturbation ip transforms into 



5(f) 

ip ->■ V'com = V + Hot = ip + U-j- 



The quantity 



11 = 4)+%- 



8<t> 



is the co-moving curvature perturbation. This quantity is gauge invariant by construction and is related to 
the gauge-dependent curvature perturbation i/>ona generic slicing to the inflaton perturbation 5(f) in that 
gauge. By construction, the meaning of 1Z is that it represents the gravitational potential on co-moving 
hypersurfaces where 5<p = or the inflaton fluctuation hypersurfaces where ip = 0: 



V>=o 



ft = V^ =0 = H 

The power spectrum of the curvature perturbation may then be easily computed 

5<f>k 



H 



(76) 



We may now compute the power spectrum of the co-moving curvature perturbation on superhorizon 
scales 



Vn{k) = 



H\ 2 



k Y 1 *- 1 



2m P i 2 e V2ti7 \aH 



= A 2 



.OH) 



where we have defined the spectral index n-jz of the co-moving curvature perturbation as 



dlnVn 

n-jz-l= ^— ^ = 3 - 2v = 2rj - 6e. 
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We conclude that inflation is responsible for the generation of adiabatic/curvature perturbations with 
an almost scale-independent spectrum. To compute the spectral index of the spectrum Vn(k) we have 
proceeded as follows: first solve the equation for the perturbation <50k in a de Sitter stage, with H = 
constant (e = rj = 0), whose solution is Eq. ( [56] ) and then taking into account the time-evolution of 
the Hubble rate and of <fi introducing the subscript in Hk and 4>k- The time variation of the latter is 
determined by 



din <f>k 
din k 




1 = -8 = e - 77. (77) 



Correspondingly, ^ is the value of the time derivative of the inflaton field when a given wavelength 
~ k^ 1 crosses the horizon (from that point on the fluctuations remains frozen in). The curvature pertur- 
bation in such an approach would read 



-7-, H k 1 (Hi 



<fik 2vr V 4>k 



Correspondingly 



dlnVn dlnHf dhxcfo 



72 dink dink dink 

During inflation the curvature perturbation is generated on superhorizon scales with a spectrum 
which is nearly scale invariant fl3) , that is, is nearly independent of the wavelength A = ir/k: the 
amplitude of the fluctuation on superhorizon scales does not (almost) depend upon the time at which the 
fluctuation crosses the horizon and becomes frozen in. The small tilt of the power spectrum arises from 
the fact that the inflaton field is massive, giving rise to a non-vanishing r\ and because during inflation 
the Hubble rate is not exactly constant, but nearly constant, where 'nearly' is quantified by the slow -roll 
parameters e. 

Comment : From what we have found so far, we may conclude that on superhorizon scales the 
co-moving curvature perturbation 1Z and the uniform-density gauge curvature £ satisfy on superhorizon 
scales the relation 

7t k ~ 0. 

7.4 Gravitational waves 

Quantum fluctuations in the gravitational fields are generated in a similar fashion to that of the scalar 
perturbations discussed so far. A gravitational wave may be viewed as a ripple of space-time in the FRW 
background metric ([T]) and in general the linear tensor perturbations may be written as 

9/j.p = a 2 (r) [-dr 2 + (<% + hij) dx l dx j ] , 

where \hij\ <C 1. The tensor hij has six degrees of freedom, but, as we studied in Subsection 7.1, the 
tensor perturbations are traceless, S^hij = 0, and transverse d l hij = (i = 1,2,3). With these four 
constraints, there remain two physical degrees of freedom, or polarizations, which are usually indicated 
A = +, x. More precisely, we can write 

= h + e± + h x e£, 

where e + and e x are the polarization tensors which have the following properties 
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ey(-k, A) = e*.(k, A), £ e*-(k, A)e«(k, A) = 4. 

A 

Notice also that the tensors hij are gauge-invariant and therefore represent physical degrees of freedom. 
If the stress-energy momentum tensor is diagonal, as the one provided by the inflaton potential 
= dfj,4>d v (f) — 9fivC, the tensor modes do not have any source in their equation and their action can 
be written as 



-9 ijdahij d a hij, 



that is the action of four independent massless scalar fields. The gauge-invariant tensor amplitude 



t> k = ampi— /i k , 



satisfies therefore the equation 



v"+ k 



v k = 0, 

which is the equation of motion of a massless scalar field in a quasi-de Sitter epoch. We can therefore 



make use of the results present in Subsection 6.5 and Eq. (63 1 to conclude that on superhorizon scales 
the tensor modes scale like 



where 



(-)(- 



vt — £■ 

2 

Since fluctuations are (nearly) frozen in on superhorizon scales, a way of characterizing the tensor per- 
turbations is to compute the spectrum on scales larger than the horizon 



A 

This gives the power spectrum on superhorizon scales 



1,3 l.3 

ft* X — > i , 1 9 (V i ,o 

^£|/ lk | 2 = 4x2^M 2 

A 



(78) 



Vr(k) 



f —\ 



mpi 2 \2tt J \aH 



k 



Tlx 



A 



T 



aH 



Tlx 



where we have defined the spectral index nx of the tensor perturbations as 



(RuVt 
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The tensor perturbation is almost scale-invariant. Notice that the amplitude of the tensor modes depends 
only on the value of the Hubble rate during inflation. This amounts to saying that it depends only on the 
energy scale V 1 ^ associated to the inflaton potential. A detection of gravitational waves from inflation 
will therefore be a direct measurement of the energy scale associated to inflation. 



7.5 The consistency relation 

The results obtained so far for the scalar and tensor perturbations allow one to predict a consistency 
relation which holds for the models of inflation addressed in these lectures, i.e., the models of inflation 
driven by one-single field cf). We define the tensor-to-scalar amplitude ratio to be 

^qAj, TOO^ ( 2 tt mpi ) 

25 V ' y 2 7r mpi J 

This means that 



tit 

r = "T 



One-single models of inflation predict that during inflation driven by a single scalar field, the ratio be- 
tween the amplitude of the tensor modes and that of the curvature perturbations is equal to minus one-half 
of the tilt of the spectrum of tensor modes. If this relation turns out to be falsified by the future measure- 
ments of the CMB anisotropics, this does not mean that inflation is wrong, but only that inflation has not 
been driven by only one field. 



7.6 From the inflationary seeds to the matter power spectrum 

As the curvature perturbations enter the causal horizon during radiation- or matter-domination, they 
create density fluctuations <5pk via gravitational attractions of the potential wells. The density contrast 

P 



= can ^ e deduced from the Poisson equation 



fc 2 $k , n r , n $Pk- 3 2 5pk 

— — = -4ttG Opk = -47rG p = -H 
n pip 

where p is the background average energy density. This means that 



2 ( k N 



3 \aH 



k- 



From this expression we can compute the power spectrum of matter density perturbations induced by 
inflation when they re-enter the horizon during matter-domination: 

from which we deduce that matter perturbations scale linearly with the wave-number and have a scalar 
tilt 



n = n-jz = 1 + 2rj — 6e. 
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The primordial spectrum Vs p is of course reprocessed by gravitational instabilities after the uni- 
verse becomes matter-dominated. Indeed, as we have seen in Section 6, perturbations evolve after en- 
tering the horizon and the power spectrum will not remain constant. To see how the density contrast 
is reprocessed we have first to analyse how it evolves on superhorizon scales before horizon-crossing. 
We use the following trick. Consider a flat universe with average energy density p. The corresponding 
Hubble rate is 

3 H 

A small positive fluctuation dp will cause the universe to be closed: 



H 2 = — (p + dp) 



2 ' 



Substracting the two equations we find 



5p 3 k J a 2 RD 



p 8irGa 2 p \ a MD 

Notice that ~ 5pa 2 /k 2 ~ (dp/p)pa 2 /k 2 = constant for both RD and MD which confirms our 
previous findings. 

When the matter densities enter the horizon, they do not increase appreciably before matter- 
domination because before matter-domination pressure is too large and does not allow the matter in- 
homogeneities to grow. On the other hand, the suppression of growth due to radiation is restricted to 
scales smaller than the horizon, while large-scale perturbations remain unaffected. This is why the hori- 
zon size at equality sets an important scale for structure growth: 



k EQ = H- 1 (oeq) - 0.08 hMpc- 1 . 

Therefore, perturbations with k 3> k E Q are perturbations which have entered the horizon before matter- 
domination and have remained nearly constant till equality. This means that they are suppressed with 
respect to those perturbations having k <C &eq by a factor (oent/oeq) 2 = (&eqA) 2 - If we define the 
transfer function T(k) by the relation 7£fi na i = T{k) "^initial we find therefore that roughly speaking 



T(k) = 

The corresponding power spectrum will be 

V 5p {k) r* 



1 k <C k E Q, 
(k EQ /k) 2 k > k EQ . 



Of course, a more careful computation needs to include many other effects such as neutrino free- 
streaming, photon diffusion and the diffusion of baryons along with photons. It is encouraging, however, 
that this rough estimate turns out to be confirmed by present data on large-scale structures [4]. 

The next step would be to investigate how the primordial perturbations generated by inflation flow 
into the CMB to produce their anisotropics. 



8 From inflation to large-angle CMB anisotropy 

As we have already mentioned, the high temperature of the early universe maintained a low equilibrium 
fraction of neutral atoms, and a correspondingly high number density of free electrons. Coulomb scat- 
tering between the ions and electrons kept them in local kinetic equilibrium, and Thomson scattering of 
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photons tended to maintain the isotropy of the CMB in the baryon rest frame. As the universe expanded 
and cooled, the dominant element hydrogen started to recombine when the temperature fell below ~ 
4000 K. This is a factor of 40 lower than might be anticipated from the 13.6 eV ionization potential of 
hydrogen, and is due to the large ratio of the number of photons to baryons. Through recombination, the 
mean-free path for Thomson scattering grew to the horizon size and CMB photons "decoupled" from 
matter. More precisely, the probability density that photons last scattered at some time defines the vis- 
ibility function. This is now known to peak 380 kyr after the Big Bang with a width ~ 120 kyr. Since 
then, CMB photons have propagated relatively unimpeded for 13.7 Gyr, covering a co-moving distance 
~ 14.1 Gpc. The distribution of their energies carries the imprint of fluctuations in the radiation temper- 
ature, the gravitational potentials, and the peculiar velocity of the radiation where they last scattered, as 
the temperature anisotropics that we observe today. 

Temperature fluctuations in the CMB arise due to various distinct physical effects: our peculiar 
velocity with respect to the cosmic rest frame; fluctuations in the gravitational potential on the last scat- 
tering surface; fluctuations intrinsic to the radiation field itself on the last scatteringsurface; the peculiar 
velocity of the last scatteringsurface and damping of anisotropics if the universe should be re-ionized 
after decoupling. The first effect gives rise to the dipole anisotropy. Finally, there is the contribution 
from the evolution of the anisotropics from the last scattering surface till today (which we shall neglect 
from now on). 

The second effect, known as the Sachs-Wolfe effect is the dominat contribution to the anisotropy 
on large-angular scales, » #hor ~ 1°- The last three effects provide the dominant contributions to 
the anisotropy on small-angular scales, <C 1°. 

8.1 Sachs- Wolfe plateau 

We consider first the temperature fluctuations on large-angular scales that arise due to the Sachs-Wolfe 
effect. These anisotropics probe length scales that were superhorizon-sized at photon decoupling and 
therefore insensitive to microphysical processes. On the contrary, they provide a probe of the original 
spectrum of primeval fluctuations produced during inflation. 

To proceed, we consider the CMB anisotropy measured at positions other than our own and at 
earlier times. This is called the brightness function @(t, x, n) = 5T(t, x, n)/T(i). The photons with 
momentum p in a given range d 3 p have intensity I proportional to T 4 (t, x, n) and therefore 51 /I = 40. 
The brightness function depends upon the direction n of the photon momentum or, equivalently, on the 
direction of observation e = — n. Because the CMB travels freely from the last-scattering, we can write 



where xls = — a?LS n is the point of the origin of the photon coming from the direction e. The co-moving 
distance of the last scatteringdistance is xls = 2/i?o- The first term corresponds to the anisotropy 
already present at last scattering and the second term is the additional anisotropy acquired during the 
travel towards us, equal to minus the fractional pertubation in the redshift of the radiation. Notice that 
the separation between each term depends on the slicing, but the sum does not. 

Consider the redshift perturbation on co-moving slicing. We imagine the universe populated by co- 
moving observers along the line of sight. The relative velocity of adjacent co-moving observers is equal 
to their distance times the velocity gradient measured along n of the photon. In the unperturbed universe, 
we have u = Hr, leading to the velocity gradient Uij = dui/drj = u,- L j = H(t)5ij with zero vorticity 
and shear. Including a peculiar velocity field as perturbation, u = Hr + v and = H(t)Sij + \ 
The corresponding Doppler shift is 




T 



da 



a 
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The perturbed FRW equation is 



6H= ^V-v, 

while 

(S p y = -3pSH - 3HSp. 

Instead of 5p, let us work with the density contrast 5 = 5p/p. Remembering that p ~ a -3 , we find that 
6 = —35H, which gives 

V-v = -4- 

From the Euler equation ii = —p~ 1 X7p — V3>, we deduce v + Hv = — V3> — p~ 1 Vp. Therefore, 
for a ~ t 2 / 3 and negligible pressure gradient, since the gravitational potential is constant, we find 

v = -£V$ 

leading to 

The photon trajectory is adx/dt = n. Using a ~ i 2 / 3 gives 

,(,)= t d JL= 3 (<«L-t 
Jt a \t a 



Integrating by parts Eq. (79), we finally find 



ST \ 1 

— J = - [*(xls) - HO)} + e • [v(0, t ) - v(x LS , t LS )] • 

The potential at our position contributes only to the unobservable monopole and can be dropped. On 
scales outside the horizon, v = — iV$ ~ 0. The remaining term is the Sachs-Wolfe effect 

5T(e) 1 , . 1, . 
= -$(x LS ) = -^(x LS ). 

This relation has been obtained as follows. The co-moving curvature perturbation is given during the 
radiation phase by 1Z = ip+H5p/p = if)— l/3Sp 7 /p T Einstein equations set ip = $and5p 7 /p 7 = — 2<3? 
on super-horizon scales. Therefore 1Z = 5/3<I> beyond the horizon. 

At large angular scales, the theory of cosmological perturbations predicts a remarkably simple 
formula relating the CMB anisotropy to the curvature perturbation generated during inflation. 

We have seen previously that the temperature anisotropy is commonly expanded in spherical 
harmonics ^(xo,ro,n) = Ylim a e,m(xo)Yi m (n) , where xo and to are our position and the preset 
time, respectively, n is the direction of observation, £'s are the different multipoles, and (ag m a^, m ,) = 
8e,e'3m,m'Ce, where the deltas are due to the fact that the process that created the anisotropy is statisti- 
cally isotropic. The C/s are the so-called CMB power spectrum. For homogeneity and isotropy, the C/s 
are neither a function of xq, nor of m. The two-point-correlation function is related to the C/'s according 
toEq. d23l. 
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For adiabatic perturbations we have seen that on large scales, larger than the horizon on the last 
scatteringsurface (corresponding to angles larger than #hor ~ 1°) ST/T = |$(xls)- In Fourier trans- 
form 

5T (Kr ,n) = 1 ik . n(T0 _ TLs) _ 



Using the decomposition 

oo 

exp(i k • n(r - r LS )) = ^{21 + l)i e j e (k(T - T LS ))P e (k • n) 

£=0 

where ji is the spherical Bessel function of order t and substituting, we get 

/ <$r(zo,TQ,n) 6T(xq,tq,ii') 
\ T T 
1 

~ V 



(80) 



(81) 



(82) 



^^^O'o-To.n) ST{ en. r n .n') 



(27T) 



T T 
/<5T(k,r ,n) /<5T(k, r , n') \ * 



if 1 00 

^ ' ^ ' e,e>=o 

jHHro ~ r LS ))Pe(k • n)i>(k' ■ n')) 



)) 
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Inserting P^(k-n) 



4-7 T 

2£H 



T E m ^m(^)^"i( n ) an( * analogously for Pe(k' ■ n'), integrating over the direc- 



tions dn k generates 8 u ,5 mm , J2 m Y e m ( n ) Y im(n')- Using as well J2 m Y tm( n ) Y tm( n ') = H£r p e( nmn ')> 
we get 



^5T(x ,T ,n) ST(x ,T ,n') 



T 
21 + 1 

4ir 



(84) 



P,(n.n')' | f (^| 2 )fc 3 j|(A:(ro-r L s)). 



Comparing this expression with that for the Ce, we get the expression for the C^ D , where the suffix 
"AD" stands for adiabatic: 



C, 



AD 



k \9 



d>| 2 )fc 3 j|(fe(r -r L s)) 



(85) 



which is valid for 2 < £ <C (tq — tls)/tls ~ 100. 

If we generically indicate by (|<3?k| 2 )fc 3 = A 2 (kro) n ~ l , we can perform the integration and get 



e{e + i)cf D 



2ir 



n— ) - 



2 1 + ; r(^)r(£+^) 



a 2 (H.\ n - 1 



(86) 



For n ~ 1 and £ S> 1, we can approximate this expression to 

l(l+l)Cf D _ 

_ _ 



2tt 



(87) 



This result shows that inflation predicts a very flat spectrum for low i. Furthermore, since inflation 
predicts $k = f^io we find that 

tt£(£+1)C{ 
WMAP5 data imply that l{ - l+l ^ u _ 1Q -io Q , 
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(88) 
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1/4 



6.7 x 10 16 GeV 



8.2 Acoustic peaks 

To be able to calculate the power spectrum of the anisotropics even on angular scales larger than 1°, we 
need to consider the evolution of the photon anistropies. As we already mentioned, before recombination 
Thomson scattering was very efficient. As a result it is a good approximation to treat photons and baryons 
as a single fluid. This treatment is called the tight-coupling approximation and will allow us to evolve 
the perturbations until recombination. 

The equation for the photon density perturbations for one Fourier mode of wave-number k is that 
of a forced and damped harmonic oscillator 



e r t 2 2 r 

7 + — | — + * c s"7 

R ■ 1 



(1 + R) 



4 

«5 7 = --kvj + 4ip. (89) 

The photon-baryon fluid can sustain acoustic oscillations. The inertia is provided by the baryons, while 
the pressure is provided by the photons. The sound speed is c 2 s = 1/3(1 + R), with R = 3pfj/4p 7 = 
31.5 (J7f,/i 2 )(T/2.7)~ 4 [(l + z)/10 3 ] -1 . As the baryon fraction goes down, the sound speed approaches 
c 2 s — > 1/3. The third equation above is the continuity equation. 



As a toy problem, we shall solve Eq. (89 1 under some simplifying assumptions. If we con- 
sider a matter-dominated universe, the driving force becomes a constant, F = — 4/3& 2 <I>, because the 
gravitational potential remains constant in time. We neglect anisotropic stresses so that tp = $, and, 



furthermore, we neglect the time dependence of R. Equation ( 89 ) becomes that of a harmonic oscillator 
that can be trivially solved. This is a very simplified picture, but it captures most of the relevant physics 
we want to discuss. 

To obtain the final solution we need again to specify the initial conditions, we shall restrict our- 
selves to adiabatic initial conditions, the most natural outcome of inflation. In our context this means 
that initially $ = ip = <3?o, ^7 = — 8/3<I>o, and v 1 = 0. We have denoted $0 the initial amplitude of the 
potential fluctuations. We shall take $0 to be a Gaussian random variable with power spectrum P$ . 

We have made enough approximations that the evaluation of the sources in the integral solution 
has become trivial. The solution for the density and velocity of the photon fluid at recombination is 

^ + $J|ls = Y(l + 312)cos(fcc a 7Ls)-*ol2, 

%1-tls = -$o(l + 3-R)c< j sin(/cc s T LS ). (90) 



Equation (90 1 is the solution for a single Fourier mode. All quantities have an additional spatial de- 
pendence (e*^ x ), which we have not included in order to make the notation more compact. With that 
additional term the solution we have is 

'^(n) = e ikDLSCOs9 S 



T 



_ (l + 3i?) r 3R 
S = $0- — t, -[cos(kc s T LS ) 



1 + 3R)' 
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-iy 1 + R cos9sm(kc s T L s)}, (91) 

where we have neglected the $ on the left-hand side because it is a constant. We have introduced 
cos 6, the cosine of the angle between the direction of observation and the wavevector k; for example, 
k • x = /cDls cos 9 . The term proportional to cos 8 is the Doppler contribution. 

Once the temperature perturbation produced by one Fourier mode has been calculated, we need 
to expand it into spherical harmonics. The power spectrum of temperature anisotropics is expressed in 
terms of the a; m coefficients as Cti = Ylm \ a em\ 2 - The contribution to Cti from each Fourier mode is 
weighted by the amplitude of primordial fluctuations in this mode, characterized by the power spectrum 
of $0 = 3/572-, -P$ = Ak~ 3 as dictated by inflation. In practice, fluctuations on angular scale £ receive 
most of their contributions from wavevectors around k* = £/Dls, so roughly the amplitude of the power 
spectrum at multipole £ is given by the value of the sources in Eq. ( J90| ) at k*. 

After summing the contributions from all modes, the power spectrum is roughly given by 



£{£ + l)C Tl « A{[ {l + * R) cos(hc s T LS ) - R} 2 + 

(1+ 3 3fi)2 C ^in 2 (fc tCs r LS )}. (92) 

Equation ( |92] > can be used to understand the basic features in the CMB power spectra. The baryon drag 
on the photon-baryon fluid reduces its sound speed below 1/3 and makes the monopole contribution 
dominant (the one proportional to cos(fc*c s TLs)- Thus, the Cti spectrum peaks where the monopole 
term peaks, A;*c s tls = tt, 27r, 3tt, • • • , which correspond to £ pea k = nirD^s/csTLS- 

It is very important to understand the origin of the acoustic peaks. In this model the universe 
is filled with standing waves; all modes of wave-number k are in phase, which leads to the oscillatory 



terms. The sine and cosine in Eq. (92 1 originate in the time dependence of the modes. Each mode £ 
receives contributions preferentially from Fourier modes of a particular wavelength k* (but pointing in 
all directions), so to obtain peaks in Q, it is crucial that all modes of a given k be in phase. If this is not 
the case, the features in the Cti spectra will be blurred and can even disappear. This is what happens 
when one considers the spectra produced by topological defects. The phase coherence of all modes of 
a given wave-number can be traced to the fact that perturbations were produced very early on and had 
wavelengths larger than the horizon during many expansion times. 

There are additional physical effects we have neglected. The universe was radiation dominated 
early on, and modes of wavelength smaller and bigger than the horizon at matter-radiation equality 
behave differently. During the radiation era the perturbations in the photon-baryon fluid are the main 
source for the gravitational potentials which decay once a mode enters into the horizon. The gravitational 
potential decay acts as a driving force for the oscillator in Eq. ([89]), so a feedback loop is established. As 
a result, the acoustic oscillations for modes that entered the horizon before matter-radiation equality have 
a higher amplitude. In the Cti spectrum the separation between modes that experience this feedback and 
those that do not occurs at I ~ -Dls/tls- Larger £ values receive their contributions from modes that 
entered the horizon before matter-radiation equality. Finally, when a mode is inside the horizon during 
the radiation era the gravitational potentials decay. 

There is a competing effect, Silk damping, that reduces the amplitude of the large-/ modes. The 
photon-baryon fluid is not a perfect fluid. Photons have a finite mean free path and thus can random-walk 
away from the peaks and valleys of the standing waves. Thus perturbations of wavelength comparable 
to or smaller than the distance the photons can random-walk get damped. This effect can be modelled 
by multiplying Eq. 0by exp(— k 2 /k 2 ), with k s 1 oc T^(Qi,h 2 ) x l 2 . Silk damping is important for 
multipoles of order ^siik ~ &s-Dls- Finally, the last scatteringsurface has a finite width. Perturbations 
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with wavelength comparable to this width get smeared out due to cancellations along the line of sight. 
This effect introduces an additional damping with a characteristic scale k^ 1 oc <5tls- 

The location of the first peak is by itself a measurement of the geometry of the universe. In fact, 
photons propagating on geodesies from the last scattering surface to us feel the spatial geometry, whose 
properties we learned are dictated by £Iq. In fact, the location of the first peak is given by l\ ~ 220/-v/Oo- 
WMAP5 gives = l-OOl^ 03- This tells us that the spatial (local) geometry of the universe is flat. This 
is precisely what inflation predicts. 

8.3 The polarization of the CMB anisotropics 

The anisotropy field is characterized by a 2 x 2 intensity tensor Iij. For convenience, we normalize this 
tensor so that it represents the fluctuations in units of the mean intensity (ijj = 51 /Iq). The intensity 
tensor is a function of direction on the sky, n, and two directions perpendicular to n that are used to define 
its components (ei,e2). The Stokes parameters Q and U are defined as Q = (In— /22)/4and U = Ivil% 
while the temperature anisotropy is given by T = (In + l22)/4 (the factor of 4 relates fluctuations in the 
intensity with those in the temperature, / oc T 4 ). When representing polarization using "rods" in a map, 
the magnitude is given by P = \J Q 2 + U 2 , and the orientation makes an angle a = \ arctan(C//Q) 
with ei . In principle the fourth Stokes parameter V that describes circular polarization is needed, but 
we ignore it because it cannot be generated through Thomson scattering, so the CMB is not expected 
to be circularly polarized. While the temperature is invariant under a right-handed rotation in the plane 
perpendicular to direction n, Q and U transform under rotation by an angle ip as 

(Q±iU)'(n) = e^(Q±iU)(n), (93) 

where = cos tp ei + sin ip e 2 and e' 2 = — sin ip ei + cos tp e 2 . The quantities Q ±iU are said to be 
spin 2. 

We already mentioned that the statistical properties of the radiation field are usually described 
in terms of the spherical harmonic decomposition of the maps. This basis, basically the Fourier basis, 
is very natural because the statistical properties of anisotropics are rotationally invariant. The standard 
spherical harmonics are not the appropriate basis for Q ± iU because they are spin-2 variables, but 
generalizations (called ± 2 Yi m ) exist. We can expand 

(Q±iU)(n) = ^2a ±2/m±2 Y im (n). (94) 

Here Q and U are defined at each direction n with respect to the spherical coordinate system (e$, e^,). 
To ensure that Q and U are real, the expansion coefficients must satisfy a*_ 2 im = a^ t-m. The equiv- 
alent relation for the temperature coefficients is a* T im = ax/- m - Instead of a±2,e m , it is convenient to 
introduce their linear combinations a E/m = -(a 2 / m + a-2,e m )/2 and a B/m = i(a 2 ,e m ~ a_2/m)/2. 
We define two quantities in real space, E(n) = J2e, m a £/m^m(n) and B(n) = J2e, m a B ,e m Ye m (n-)- 
Here E and B completely specify the linear polarization field. 

The temperature is a scalar quantity under a rotation of the coordinate system, T'(n' = IZn) = 
T(n), where 1Z is the rotation matrix. We denote with a prime the quantities in the transformed coor- 
dinate system. While Q ± iU are spin 2, E(n) and B(n) are invariant under rotations. Under parity, 
however, E and B behave differently, E remains unchanged, while B changes sign. 

To characterize the statistics of the CMB perturbations, only four power spectra are needed, those 
for T, E, B and the cross correlation between T and E. The cross correlation between B and E or B 
and T vanishes if there are no parity-violating interactions because B has the opposite parity to T or 
E. The power spectra are defined as the rotationally invariant quantities Cxe = 2E+1 Em(4 e m a T,em) , 

C E e = 2l+l^rri( a *E,em a E,£rn}> C B i = 5^4 Em( a B/m fl B/m)' ^ C Cl = 2l+l Y. m ( a T,e m a E,em) ■ 

The brackets (• • • ) denote ensemble averages. 
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Fig. 4: Examples of E- and £>-mode patterns of polarization 
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Fig. 5: Thomson scattering of radiation where quadrupole anisotropy generates linear polarization 

Polarization is generated by Thomson scattering between photons and electrons, which means that 
polarization cannot be generated after recombination (except for re-ionization, which we shall discuss 
later). But Thomson scattering is not enough. The radiation incident on the electrons must also be 
anisotropic. In fact, its intensity needs to have a quadrupole moment. This requirement of having both 
Thomson scattering and anisotropies is what makes polarization relatively small. After recombination, 
anisotropics grow by free streaming, but there is no scattering to generate polarization. Before recom- 
bination there were so many scatterings that they erased any anisotropy present in the photon-baryon 
fluid. 

In the context of anisotropies induced by density perturbations, velocity gradients in the photon- 
baryon fluid are responsible for the quadrupole that generates polarization. Let us consider a scattering 
occurring at position xq: the scattered photons came from a distance of order the mean free path (Ay) 
away from this point. If we are considering photons traveling in direction h, they roughly come from x = 
xq + Xrn. The photon-baryon fluid at that point was moving at velocity v(x) « v(xq) + \Tnid{v(xo). 
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Due to the Doppler effect the temperature seen by the scatterer at xq is 5T(xQ,h) = h-[v(x) —v(xq)] « 
\TnihjdiVjixo), which is quadratic in h (i.e., it has a quadrupole). Velocity gradients in the photon- 
baryon fluid lead to a quadrupole component of the intensity distribution, which, through Thomson 
scattering, is converted into polarization. 

The polarization of the scattered radiation field, expressed in terms of the Stokes parameters Q 
and U, is given by (Q + ill) oc <jt f d£l'(m • h') 2 T(h') oc Xpia^mP diVj\ TLS , where ot is the Thomson 
scattering cross-section and we have written the scattering matrix as P(m, n') = — 3/4crr(m • n') 2 , with 
m = ei + 1&2 ■ In the last step, we integrated over all directions of the incident photons n'. As photons 
decouple from the baryons, their mean free path grows very rapidly, so a more careful analysis is needed 
to obtain the final polarization: 

(Q + iU)(n) rj eST hS m i m j diV j \ TLS , (95) 

where <5tls is the width of the last scattering surface and gives a measure of the distance that photons 
travel between their last two scatterings, and e is a numerical constant that depends on the shape of the 



visibility function. The appearance of m l m J in Eq. (95) ensures that (Q + ill) transforms correctly 
under rotations of (ei, 62). 

If we evaluate Eq. ( |93| ) for each Fourier mode and combine them to obtain the total power, we get 
the equivalent of Eq. d92b, 



£{i + l)C Ei « Ae 2 (l + 3R) 2 (K5t ls ) 2 sin 2 (Awls), (96) 
where we are assuming n = 1 and that £ is large enough that factors like (I + 2)\/(£ — 2)1 « £ 4 . 



The extra fc* in Eq. ( |96| ) originates in the gradient in Eq. (95 1. The large-angular scale polarization 
is greatly suppressed by the Mtls factor. Correlations over large angles can only be created by the 
long-wavelength perturbations, but these cannot produce a large polarization signal because of the tight 
coupling between photons and electrons prior to recombination. Multiple scatterings make the plasma 
very homogeneous; only wavelengths that are small enough to produce anisotropics over the mean free 
path of the photons will give rise to a significant quadrupole in the temperature distribution, and thus 
to polarization. Wavelengths much smaller than the mean free path decay due to photon diffusion (Silk 
damping) and so are unable to create a large quadrupole and polarization. As a result polarization peaks 
at the scale of the mean free path. 

On sub-degree angular scales, temperature, polarization, and the cross-correlation power spectra 
show acoustic oscillations. In the polarization and cross-correlation spectra the peaks are much sharper. 
The polarization is produced by velocity gradients of the photon — baryon fluid at the last scatteringsur- 
face. The temperature receives contributions from density and velocity perturbations, and the oscillations 
in each partially cancel one another, making the features in the temperature spectrum less sharp. The 
dominant contribution to the temperature comes from the oscillations in the density [Eq. ([90])], which 



are out of phase with the velocity. This explains the difference in location between the temperature 



and polarization peaks. The extra gradient in the polarization signal, Eq. (95 ), explains why its overall 
amplitude peaks at a smaller angular scale. 

Now, as photons travel in the metric perturbed by a GW [ds 2 = a 2 (r) [—dr 2 +(Sij+hfj)dx i dx^]], 
they get redshifted or blueshifted depending on their direction of propagation relative to the direction of 
propagation of the GW and the polarization of the GW. For example, for a GW travelling along the z 
axis, the frequency shift is given by 

~ = \ n l n j hl {±) = 1(1- cos 2 e)e ±i2(t ' h t exp(ik • x), (97) 
v uT 2 J 2 

where (6, (ft) describe the direction of propagation of the photon, the ± correspond to the different polar- 
izations of the GW, and h t gives the time-dependent amplitude of the GW. During the matter-dominated 
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Fig. 6: E- and B-mode power spectra for a tensor-to-scalar ratio saturating the current bounds, r = 0.3 and for 
r = 0.01. Shown are the experimental sensitivities of WMAP, Planck and two different realizations of CMBPol 
(EPIC-LC and EPIC-2m) 



era, for example, h t = 3ji(kr)/kr: time changes in the metric lead to frequency shifts (or equivalently 
shifts in the temperature of the black body spectrum). Notice that the angular dependence of this fre- 
quency shift is quadrupolar in nature. As a result, the temperature fluctuations induced by this effect 
as photons travel between successive scatterings before recombination produce a quadrupole intensity 
distribution, which, through Thomson scattering, lead to polarization. Both E and B power spectra are 
generated by GW. The current push to improve polarization measurements follows from the fact that 
density perturbations, to linear order in perturbation theory, cannot create any B-type polarization. As a 
rough rule of thumb, the amplitude of the peak in the B-mode power spectrum for GW is 



[£(£ + l)C , K /2vr] 1 / 2 = 0.024(T/ 1 / 4 /10 16 GeV) 2 /iK 



where 



V 1/4 ~ 6.7 r 1 / 4 x 10 16 GeV (98) 

is the energy scale of inflation. A future experiment like CMBPol [ 14] can probe values of r as small 
as 10 -2 , corresponding to an inflation energy scale of about 2 x xlO 16 GeV. Furthermore, using the 
consistency relation r = e valid in one-single field models of inflation, one deduces that 

^, (99) 



TOpi VlO 2 / 

meaning that a future measurement of the B-mode of CMB polarization will imply an inflaton excursus 
of Planckian values. Therefore, A future measurement of the B-mode polarization of the CMB will 
allow a determination of the value of the energy scale of inflation. This explains the utility of CMB 
polarization measurements as probes of the physics of inflation. A detection of primordial B-mode 
polarization would also demonstrate that inflation occurred at a very high energy scale, and that the 
inflaton traversed a super-Planckian distance in field space. 



9 The dark puzzles 

Having explored the physics of the primordial epochs of the evolution of the universe, such as inflation, 
and its impact on the present-day observables, we now devote the remaining space to a short discussion 
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of the dark puzzles of the present-day universe: the dark energy and the dark matter puzzles. 



9.1 A present-day accelerating universe 

In 1998 the accelerated expansion of the universe was pointed out by two groups from the observations 
of Type la Supernova (SN la) p3|[16| . Let us see how this came about. 

An important concept related to observational tools in an expanding background is associated with 
the definition of a distance. A way of defining a distance is through the luminosity of a stellar object. 
The distance <1l known as the luminosity distance, plays a very important role in astronomy including in 
supernovae observations. It proves to be convenient to write the metric as 

ds 2 = -dt 2 + a 2 {t) [dr 2 + f 2 K {r){d6 2 + sin 2 6<# 2 )] , (100) 

where 

sinr , K = +1 , 
I Kir) - { r, K = 0, (101) 

sinhr , K = — 1 . 

In Minkowski space time the absolute luminosity L s of the source and the energy flux J 7 at a distance d is 
related through F = L s / (And 2 ). By generalizing this to an expanding universe, the luminosity distance, 
di, is defined as 

d\ = -A- . (102) 

Let us consider an object with absolute luminosity L s located at a co-moving coordinate distance r from 
an observer at r = 0. The energy of light emitted from the object with time interval At e is denoted as 
AE e , whereas the energy which reaches at the sphere with radius r is written as AE r . We note that AE e 
and AE r are proportional to the frequencies of light at r and r = 0, respectively, i.e., AE e oc v e and 
AE r oc v r . The luminosities L r and L e are given by 

AE e AE e 
L r = —^, L e = — t. (103) 
At e At e 

The speed of light is given by c = v e \ e = v r \ r , where A e and A r are the wavelengths at r and r = 0. 
Then, we find 

A r v e At r AE e 
A e v r At e AE r 
where we have also used v r At r = v e At e - Combining Eq. ( 103 1 with Eq. ( |104 i, we obtain 

L e = L r {l + z) 2 . (105) 

The light travelling along the r direction satisfies the geodesic equation ds 2 = —dt 2 + a 2 (t)dr 2 = 0. 
We then obtain 

r f u dt 



J 7 = i—, — ; , ^ 2 • (107) 



From the metric ( 100 1 we find that the area of the sphere at t = t r is given by S = 4ir(a r fK(r)) 2 . Hence 
the observed energy flux is 

L 

An(a r f K {r)) 

Substituting Eqs. (106 1 and ( 107 > for Eq. ( |102| ), we obtain the luminosity distance in an expanding 
universe: 

d L = a r f K (r)(l + z). (108) 
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In the flat FRW background with /ftr(r) = r we find 



l + z 



dz 



i H 



H J J - H(z'Y 
Then the Hubble rate H(z) can be expressed in terms of cIl(z): 

H(z) f d {dL{z) 



dz \ 1 + z 



(109) 



(110) 



If we measure the luminosity distance observationally, we can determine the expansion rate of the uni- 
verse. 

The energy density p on the right-hand side of Einstein equations includes all components present 
in the universe, namely, non-relativistic particles, relativistic particles, cosmological constant and so on 



P = (a/°o) 



-3(1+Wi) 



Y^ P f\i + zf^. 



(in) 



Here W{ and pf* correspond to the equation of state and the present energy density of each component, 

(112) 



respectively. The Hubble parameter takes the convenient form 

H 2 = HlY J ^f\ 1 + z ? {1+Wl \ 



where tlf' = 8ttGp[°' / (SHq) = pf^/p^ is the density parameter for an individual component at the 
present epoch. Hence the luminosity distance in a flat geometry is given by 



[l + z) 



H, 



o Jo 



dz' 



(113) 



The direct evidence for the current acceleration of the universe is related to the observation of luminosity 
distances of high redshift supernovae | [T5} 16 1. The apparent magnitude m of the source with an absolute 



magnitude M is related to the luminosity distance d^ via the relation 1 17 1 



m 



M = 5 log 



10 



Mpc 



+ 25. 



(114) 



This comes from taking the logarithm of Eq. ( 102| ) by noting that m and M are related to the logarithms 
of F and L s , respectively. The numerical factors arise because of conventional definitions of m and M 
in astronomy. Type la supernovae (SN la) can be observed when white dwarf stars exceed the mass of 
the Chandrasekhar limit and explode. The belief is that SN la are formed in the same way irrespective of 
where they are in the universe, which means that they have a common absolute magnitude M independent 
of the redshift z. Thus they can be treated as an ideal standard candle. We can measure the apparent 
magnitude m and the redshift z observationally, which of course depends upon the objects we observe. 

In order to get a feeling of the phenomenon let us consider two supernovae 1992P at low-redshift 
z = 0.026 with m = 16.08 and 1997ap at high-redshift z = 0.83 with m = 24.32 (BJ. As we 
have already mentioned, the luminosity distance is approximately given by di(z) ~ z/Hq for z <C 1. 
Using the apparent magnitude m = 16.08 of 19 92P at z = 0.026, we find that the absolute magnitude 
is estimated by M = -19.09 from Eq. \\U\ . Here we adopted the value Hq 1 = 2998/t^ 1 Mpc 
with h = 0.72. Then the luminosity distance of 1997ap is obtained by substituting m = 24.32 and 
M -- 



19.09 for Eq. (114 1: H d L ~ 1.16 for z 



luminosity distance in a two-component flat universe is H^di 
for d$ ~ 0.3, Ode ^ 0.7 



0.83. From Eq. ( |1 13[ ) the theoretical estimate for the 

1 and H d L ~ 1.23 



0.95 for 
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Fig. 7: The luminosity distance versus redshift for a flat cosmological model. Black points are from the "Gold'' 
data sets (18); red points are from recent data from HST 



In 2004 Riess et al. (T8 ] reported the measurement of 16 high-redshift SN la with redshift z > 1.25 
with the Hubble Space Telescope (HST). By including 170 previously known SN la data points, they 
showed that the universe exhibited a transition from deceleration to acceleration at > 99% confidence 



level. A best-fit value of was found to be O^? 



0.29 



+0.05 
-0.03 



(the error bar is la). This shows that a 



matter-dominated universe without a cosmological constant does not fit the data. 

We should emphasize that the accelerated expansion is by cosmological standards really a late- 
time phenomenon, starting at a redshift z ~ 1. From Eq. (112) the deceleration parameter, q = —ad /a 2 , 
is given by 



(z) = 3 E^f(i + ^)(i + ^) 3(1+w ' ) _ x 
2 EA (0) (i + ^) 3{1+u,i) 



For the two-component flat cosmology, the universe enters an accelerating phase (q < 0) for z < z c = 
(2$7de/^dm) 1 ^ 3 — 1- When Odm = 0.3 and Ode = 0.7, we have z c = 0.67. The problem of why 
an accelerated expansion should occur now in the long history of the universe is called the "coincidence 
problem". 

9.1.1 The origin of the acceleration 

Once the idea of the accelerating universe is accepted, the next pressing question is: Why? There are 
various explanations available that we may mention briefly. The general trend is to accept that there is a 
form of Dark Energy (DE) fluid dominating the energy density of the present day. Its pressure is P = wp 
and w needs to be smaller than —1/3 for this fluid to cause the acceleration. Having learned how to use 
scalar fields to accelerate the universe at primordial epochs, the most natural way to explain DE whould 
be to introduce a scalar field dubbed quintessence, with potential 



V(0) = V + V ((/>), 



(115) 
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Fig. 8: The dark energy (vacuum energy) and dark matter (mass density) abundances from SN, CMB, and galaxy 
clustering observations 



Now, if Vq 3> V(<fi) (at least at present epochs), the DE is in practice a Cosmological Constant (CC). 
Its value must be extremely small, ~ (Hompi) 1 / 2 ~ 10~ 3 eV. Why it is so small is a mystery 

that earned the name "the CC problem". On the other hand, if Vq <C V (</>), then the dynamics of 
the quintessence field dominates. However, another problem arises at this stage. Having learned from 
inflation that the field must be slow-rolling to cause the acceleration of the universe, we have to assume 
that (rrip^V" /V) is smaller than unity. This implies that <f> is of order of the Planck scale and that its 
mass squared is such that V" ~ Hq ~ (10~ 33 eV) 2 . The quintessence field has a Compton wavelegth 
as large as the entire observed universe. 

If the reader does not like all this fine-tuning, there are at least two other explanations for the 
acceleration of the universe. The first one goes under the name of modified gravity and is in fact rather 
intuitive. If gravity gets weaker at large distances, objects far from us may recede at a velocity larger 
than what they would do in the traditional Newtonian gravity case. For this to work, we have to suppose 
that the gravitational force has a transient at some critical (and cosmological) scale r c , from the usual 
1/r 2 to, say 1 /r 3 . How to get this transition is unfortunately beyond the scope of these lectures. Another 
alternative goes under the name of the "anthropic principle" and is based on the following point. As 
we have seen, in a static universe, overdense regions will increase their density at an exponential rate. 
In an expanding universe, however, there is a competition between the expansion and the gravitational 
collapse. More rapid expansion, as induced by DE, retards the growth of structure. General relativity 
provides the following useful relation in linear perturbation theory between the growth factor g(z) and 
the expansion history of the universe 

g + 2Hg = 4nGp m = g . (116) 
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If the universe is always matter-dominated, then g ~ a; however, in a DE dominated universe g scales 
slower than the scale factor. Now, if the CC is too large, structure does not have time to develop: the 
initial condition is Sp m /p m ~ 10~ 5 at the last scattering surface (z ~ 10 3 ) and needs to becomes order 
unity by now. Now, if we impose that structures might have been able to develop by now even in the 
presence of a CC, one obtains a reassuring bound, the CC V^ 4 must be smaller than about 10 _1 eV. In 
other words, the CC may not be far from the value we observe (if it is non-zero) because otherwise we 
would not be here to discuss about it. A great deal of observational effort of the next decades will be 



devoted to understand the cause of the acceleration of the universe [19]. Four observational techniques 
are currently receiving much attention: 1) Baryoniuc Acoustic Oscillations (BAO) are observed on large- 
scale surveys of the spatial distribution of matter. They are caused by the same oscillations that left an 
imprint in the CMB under the form of acoustic peaks. The BAO technique is sensitive to the DE through 
its effect on the angular-diameter distance vs. redshift relation and through its effect on the time evolution 
of the expansion rate; 2) Galaxy Cluster (CL) surveys measure the spatial density and distribution of 
galaxy clusters. The CL technique is sensitive to DE through its effect in the angular-diameter distance 
vs. redshift relation and through its effect on the time evolution of the expansion rate and the growth 
rate of perturbations; 3) supernovae as standard candles to determine the luminosity distance vs. redshift 
relation; 4) Weak Lensing (WL) surveys measure the distortion of background images due to the bending 
of light as it passes by galaxies or clusters of galaxies. The WL technique is sensitive to DE through its 
effect on the angular-diameter distance vs. redshift relation and the growth rate of perturbations. All 
these techniques will not only shed light on the nature of DE, but will also help us to discriminate the 
various possibilities to explain the present-day acceleration. For instance, the modified gravity scenario 
predicts a growth function which is different from the one predicted in a CC dominated universe. Future 
applications of the techniques briefly summarized above should be able to determine which scenario is 
more likely. 

9.2 Dark matter 

The evidence that 95% of the mass of galaxies and clusters is made of some unknown component of 
Dark Matter (DM) comes from (i) rotation curves (out to tens of kpc), (ii) gravitational lensing (out to 
200 kpc), and (iii) hot gas in clusters. They lead us to believe that DM makes up about 30% of the entire 



energy of the universe. A nice review about DM can be found in Ref. [20|. 

In the 1970s, Ford and Rubin discovered that rotation curves of galaxies are flat. The velocities 
of objects (stars or gas) orbiting the centres of galaxies, rather than decreasing as a function of the 
distance from the galactic centres as had been expected, remain constant out to very large radii. Similar 
observations of flat rotation curves have now been found for all galaxies studied, including our Milky 
Way. The simplest explanation is that galaxies contain far more mass than can be explained by the bright 
stellar objects residing in galactic disks. This mass provides the force to speed up the orbits. To explain 
the data, galaxies must have enormous dark haloes made of unknown matter. Indeed, more than 95% of 
the mass of galaxies consists of dark matter. The baryonic matter which accounts for the gas and disk 
cannot alone explain the galactic rotation curve. However, adding a DM halo allows a good fit to data. 

The limitations of rotation curves are that one can only look out as far as there is light or neutral 
hydrogen (21 cm), namely to distances of tens of kpc. Thus one can see the beginnings of DM haloes, 
but cannot trace where most of the DM is. The lensing experiments discussed in the next section go 
beyond these limitations. 

Einstein's theory of General Relativity predicts that mass bends, or lenses, light. This effect can 
be used to gravitationally ascertain the existence of mass even when it emits no light. Lensing measure- 
ments confirm the existence of enormous quantities of DM both in galaxies and in clusters of galaxies. 
Observations are made of distant bright objects such as galaxies or quasars. As the result of intervening 
matter, the light from these distant objects is bent towards the regions of large mass. Hence there may 
be multiple images of the distant objects, or, if these images cannot be individually resolved, the back- 
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ground object may appear brighter. Some of these images may be distorted or sheared. The Sloan Digital 
Sky Survey used weak lensing (statistical studies of lensed galaxies) to conclude that galaxies, including 
the Milky Way, are even larger and more massive than previously thought, and require even more DM 
out to great distances. Again, the predominance of DM in galaxies is observed. The key success of 
the lensing of DM to date is the evidence that DM is seen out to much larger distances than could be 
probed by rotation curves: the DM is seen in galaxies out to 200 kpc from the centres of galaxies, in 
agreement with N-body simulations. On even larger Mpc scales, there is evidence for DM in filaments 
(the cosmic web). Another piece of gravitational evidence for DM is the hot gas in clusters. The X-ray 
data indicates the presence of hot gas. The existence of this gas in the cluster can only be explained by 
a large DM component that provides the potential well to hold on to the gas. In summary, the evidence 
is overwhelming for the existence of an unknown component of DM that comprises 95% of the mass in 
galaxies and clusters. 

There is another basic reason why DM is necessary: to form structures as we observe them. Let us 
assume that the matter content of the universe is dominated by a pressureless and self-gravitating fluid. 
This approximation holds if we are dealing with the evolution of the perturbations in the DM component 
or in case we are dealing with structures whose size is much larger than the typical Jeans scale length of 
baryons. Let us also define x to be the co-moving coordinate and r = a(t)x the proper coordinate, a(t) 
being the cosmic expansion factor. Furthermore, if v = r is the physical velocity, then v = dx + u, 
where the first term describes the Hubble flow, while the second term, u = o(t)x, gives the peculiar 
velocity of a fluid element which moves in an expanding background. 

In this case the equations that regulate the Newtonian description of the evolution of density per- 
turbations are the continuity equation: 

Or 

— + V-[(l + <5)u]=0, (117) 
which gives the mass conservation, the Euler equation 



dt 



+ 2#(i)u + (u • V)u = (118) 



a 



which gives the relation between the acceleration of the fluid element and the gravitational force, and the 
Poisson equation 

V 2 (j) = 4irGpa 2 5 (119) 

which specifies the Newtonian nature of the gravitational force. In the above equations, V is the gradient 
computed with respect to the co-moving coordinate x, </>(x) describes the fluctuations of the gravitational 
potential, and H(t) = a/a is the Hubble parameter at the time t. Its time-dependence is given by 
H(t) = E(t)H , where 

e(z) = [(i + zfn m + (i + zf(i - n m - n DE ) + (i + Z )^ 1+W ^ DE ] 1/2 . (120) 

In the case of small perturbations, these equations can be linearized by neglecting all the terms 
which are of second order in the fields 5 and u. In this case, using the Euler equation to eliminate the 
term du/dt, and using the Poisson equation to eliminate V 2 4>, one ends up with 

-^ + 2H{t)—- iirGpS = . (121) 

This equation describes the Jeans instability of a pressureless fluid, with the additional "Hubble drag" 
term 2H(t)dS/dt, which describes the counter-action of the expanding background on the perturbation 
growth. Its effect is to prevent the exponential growth of the gravitational instability taking place in a 
non-expanding background. The solution of the above equation can be cast in the form: 

5(x,*) = «y+(x,t i )D+(t) + 5_(x,t i ) J D_(t), (122) 
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where D + and D_ describe the growing and decaying modes of the density perturbation, respectively. 
In the case of an Einstein-de-Sitter (EdS) universe (Q m = 1, £Ide = 0), it is H{t) = 2/(3t), so that 
D+(t) = (t/U) 2 / 3 and £>_(*) = (t/U)' 1 . The fact that £)+(*) oc a(t) for an EdS universe should not 
be surprising. Indeed, the dynamical time-scale for the collapse of a perturbation of uniform density p is 
*dyn « (Gp) -1 / 1 2 , while the expansion time-scale for the EdS model is t cxp oc (Gp)~ l l I 2 , where p is the 
mean cosmic density. Since for a linear (small) perturbation it is p ~ p, then t^yn ~ ^cxp, thus showing 
that the cosmic expansion and the perturbation evolution take place at the same pace. This argument 
also leads to understanding the behaviour for a £l m < 1 model. In this case, the expansion time scale 
becomes shorter than the above one at the redshift at which the universe recognizes that Q rn < 1. This 
happens at 1 + z ~ fij' 3 or at 1 + z ~ fi" 1 in the presence or absence of a cosmological constant 
term, respectively. Therefore, after this redshift, cosmic expansion takes place at a quicker pace than 
gravitational instability, with the result that the perturbation growth is frozen. 

The exact expression for the growing model of perturbations is given by 

r poo -1 , / 

D + {z) = -n m E(z) J ^i^dz'. (123) 

The EdS has the faster evolution, while the slowing down of the perturbation growth is more apparent for 
the open low-density model, the presence of a cosmological constant providing an intermediate degree of 
evolution. The key point is, however, that a pressureless fluid such as DM is needed for the perturbations 
to grow to give rise to collapsed objects. Baryon perturbations, being coupled to photons till the last- 
scattering epoch, feel a non-vanishing pressure and therefore they may not grow. After the last- scattering 
stage, the baryons fall into the gravitational potential generated by DM and the baryonic perturbations 
may promptly catch up with those of DM. 

9. 2. 1 Dark matter candidates 

There is a plethora of dark matter candidates. MACHOs, or Massive Compact Halo Objects, are made 
of ordinary matter in the form of faint stars or stellar remnants; they could also be primordial black holes 
or mirror matter. However, there are not enough of them to completely resolve the question. Of the 
non-baryonic candidates, the most popular are the WIMPS (Weakly Interacting Massive Particles) and 
the axions, as these particles have been proposed for other reasons in particle physics. Ordinary massive 
neutrinos are too light to be cosmologically significant, though sterile neutrinos remain a possibility. 
Other candidates include primordial black holes, non-thermal WIMPzillas, and Kaluza-Klein particles 
which arise in higher dimensional theories. 

About axions, the good news is that cosmologists do not need to "invent" new particles. Two 
candidates already exist in particle physics for other reasons: axions and WIMPs. Axions with masses 
in the range 10" ( 3 ~ 6 ) eV arise in the Peccei-Quinn solution to the strong-CP problem in the theory of 
strong interactions. 

WIMPs are also natural dark matter candidates from particle physics. These particles, if present 
in thermal abundances in the early universe, annihilate with one another so that a predictable number of 
them remain today. The relic density of these particles comes out to be the right value: 

1]dm/j 2 = (3 x 10" 26 cm 3 /s)/(du) A (124) 

where the annihilation cross-section {av)x of weak interaction strength automatically gives the right 
answer. The reason why the final abundance is inversely proportional to the annihilation cross-section 
is rather clear: the larger the annihilation cross-section, the more WIMPs annihilate and the fewer of 
them are left behind. Furthermore, annihilation is not eternal: owing to the expansion of the universe, 
annihilation stops when its rate becomes smaller than the expansion rate of the universe. When this 
happens, the abundance is said to freeze-out. 
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Fig. 9: The abundance of WIMPs of a given mass to as a function of temperature and for various annihilation 
cross-sections 



This coincidence is known as 'the WIMP miracle' and is the reason why WIMPs are taken so 
seriously as DM candidates. The best WIMP candidate is motivated by Supersymmetry (SUSY): the 
lightest neutralino in the Minimal Supersymmetric Standard Model. Supersymmetry in particle theory 
is designed to keep particle masses at the right value. As a consequence, each particle we know has a 
partner: the photino is the partner of the photon, the squark is the quark's partner, and the selectron is the 
partner of the electron. The lightest superysmmetric partner is a good dark matter candidate. 

There are several ways to search for dark WIMPs. SUSY particles may be discovered at the LHC 
as missing energy in an event. In that case one knows that the particles live long enough to escape the 
detector, but it will still be unclear whether they are long-lived enough to be the dark matter. Thus com- 
plementary astrophysical experiments are needed. In direct detection experiments, the WIMP scatters off 
a nucleus in the detector, and a number of experimental signatures of the interaction can be detected. In 
indirect detection experiments, neutrinos that arise as annihilation products of captured WIMPs exit from 
the Sun and can be detected on Earth. Another way to detect WIMPs is to look for anomalous cosmic 
rays from the Galactic Halo: WIMPs in the Halo can annihilate with one another to give rise to antipro- 
tons, positrons, or neutrinos. In addition, neutrinos, gamma rays, and radio waves may be detected as 
WIMP annihilation products from the Galactic Centre. For lack of time these issues were not discussed 
extensively in the lectures. The interested reader may find more about these issues in Ref. [20j. 



10 Conclusions 

The period when we say that cosmology is entering a golden age has already passed: cosmology is in the 
middle of its golden age. Present observational data pose various puzzles whose solutions might either be 
around the corner or decades far in the future. It will require some young and creative researcher sitting 
in this room to solve them. This is why the cosmological puzzles are dark, but the future is brighter. 
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